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ABSTRACT 


This  report  describes  the  results  of  the  continued  evaluation 
of  the  Norwegian  Long-Period  Seismic  Array  (NORSAR)  by  Texas  Instruments 
Incorporated  at  the  Seismic  Data  Analysis  Center  during  the  period  1  April  1972 

to  31  March  1973  The  NORSAR  evaluation  was  concerned  with  five  major  areas 
of  study: 

•  Noise  analysis 

•  Array  processing  performance 

•  Effectiveness  of  matched  filters 

•  Surface  wave  detection  capability 

•  Performance  of  standard  surface  wave  discriminants 

This  evaluation  used  324  Eurasian  events  and  32  noise  samples 
taken  from  the  period  of  1  January  1972  to  26  December  1972.  These  data  were 
combined  with  earlier  NORSAR  data,  when  applicable,  so  as  to  maximize  the 
data  base. 

Neither  the  Advanced  Research  Projects  Agency  nor  the  Air  Force 
Technical  Applications  Center  will  be  responsible  for  information  contained 
herein  which  has  been  supplied  by  other  organizations  or  contractors,  and  this 
document  is  subject  to  later  revision  as  may  be  necessary.  The  views  and  con¬ 
clusions  presented  are  those  of  the  authors  and  should  not  be  interpreted  as 
necessarily  representing  the  official  policies,  either  expressed  or  implied,  of 
the  Advanced  Research  Projects  Agency,  the  Air  Force  Technical  Applications 
Center,  or  the  US  Government. 
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SECTION  I 


INTRODUCTION 


This  report  presents  results  of  the  evaluation  of  the  Norwegian 
Long-Period  Seismic  Array  (NORSAR)  made  during  the  period  from  April  1972 
through  March  1973.  This  work  is  a  continuation  of  an  analysis  of  NORSAR 
reported  on  earlier  (Eyres,  Laun,  and  Swindell,  1972  ).  The  eval¬ 
uation  was  directed  toward  two  major  goals:  the  determination  of  the 

detection  and  discrimination  capabilities  of  NORSAR  and  the  investigation  of 
methods  which  sustain  or  enhance  those  capabilities.  Five  separate  studies 
were  made: 

•  Noise  analysis 

•  Array  processing  effectiveness 

•  Matched  filler  performances 

•  Signal  detection  threshold  estimation 

•  Behavior  of  standard  discriminants 

The  results  from  these  studies  are  discussed  in  Sections  III 
through  VII.  When  possible,  results  from  the  first  report  have  been  included 
with  those  obtained  this  year  so  as  to  maximize  the  data  base.  The  total  data 
base  is  described  in  Section  II.  Section  VIII  summarizes  the  results  and 
conclusions  and  suggests  areas  of  future  NORSAR  long-period  analysis. 
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SECTION  II 


DATA  BASE 


The  results  presented  in  this  report  are  based  mainly  on 
seismic  events  recorded  during  two  periods  in  1972:  1  January  through 
20  March,  and  1  June  through  31  July.  Other  events  from  1972  include  all 
reported  presumed  explosions  from  the  Sino-Soviet  area.  There  were  324 
events  successfully  processed  from  1972.  Finally,  59  events  from  1971, 
which  were  reported  on  last  year,  have  been  included. 

Table  II - 1  lists  these  383  events  and  various  parameters 
associated  with  them,  including:  name  (the  three-part  designator  consists 
of  region,  julian  date,  and  hour  of  origin  time.  A  "/"  indicates  a  1971  event; 
an  indicates  a  1972  event),  date  and  origin  time,  epicenter  location,  delta 
and  azimuth  from  NORSAR,  bodywave  magnitude  (mb),  surface  wave  magnitudes 
for  vertical  (Rayleigh)  and  transverse  (Love)  components,  converted  AR  and 
AL  values,  depth  (if  reported),  and  coded  comments.  The  comments  include 
a  letter  indicating  the  source  of  event  information,  an  :  X"  if  the  event  is  a 
presumed  explosion,  and  the  number  of  sites  used  to  form  the  array  beam. 

The  source  codes  are  "P"  for  the  National  Oceanic  and  Atmospheric  Ad¬ 
ministration  (NOAA)  Preliminary  Determination  of  Epicenters  (PDE)  bulletin, 
ML"  for  the  Montana  Large  Aperture  Seismic  Array  (LASA)  bulletin,  "N"  for 
the  NORSAR  bulletin,  and  either  an  nIn  or  a  "J"  for  the  Massachusetts  Insti¬ 
tute  of  Technology -Lincoln  Lab  bulletin  for  the  International  Seisrrological 
Month  (ISM),  which  was  from  20  February  to  19  March  1972. 
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Miscellaneous  symbols  which  appear  on  the  table  are:  an 
asterisk  following  the  m^  if  the  was  recalculated  using  only  stations  at 

teleseismic  (>15°)  distances;  a  "C"  following  an  M  if  detection  was  obtained 

s 

only  with  a  chirp  filter.  Non-detection  of  an  event  on  a  particular  component 

is  indicated  by  no  entry  in  the  M  column.  Love  wave  magnitudes  were  not 

s 

computed  for  the  1971  events. 

The  324  events  from  1972  which  were  successfully  processed 
were  obtained  from  a  list  552  events  originally  selected  for  processing.  Two 
hundred  twenty-eight  events  were  rejected  as  being  unsuitable  for  further 


processing  for  the  following  reasons: 

Lackof  data  or  excessive  time  gaps  104 

Parity  errors,  spikes,  other  tape  problems  13 

Events  subsequently  found  to  be  mislocated  3 

Interference  from  other  events  83 


Events  originally  taken  from  LASA  or  NORSAR 
bulletins  which  PDE  data  later  showed  to  be 
deep  (>50km)  .  25 

228 

Thus,  59%  of  the  events  initially  selected  for  analysis  were  ultimately 
successfully  processed. 

The  data  used  for  noise  analysis  and  array  processing  analy¬ 
sis  were  edited  quasi-periodically  at  ten  day  intervals  continuing  the  analysis  of 
1971  (Eyres  etal.,  1972).  The  intervals  were  selected  so  as  to  exclude  any  surface 

wave  arrivals  from  known  events  at  the  array.  Exceptionally  long  intervals 
of  signal-free  noise  were  edited  whenever  found.  The  data  from  at  least  one 
site  of  eachnoise  sample  were  always  plotted  to  make  sure  that  surface  waves 
from  unreported  events  were  not  present. 
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A  diagram  of  the  22 -site  NORSAR  long-period  array  is 
shown  in  Figure  II— 1 . 

In  contrast  to  the  various  data  quality  problems  of  last 
year,  the  data  from  the  first  seven  months  of  1972  were  of  generally  high 
quality.  This  year  an  average  of  16.  5  sites  were  used  for  beam  forming  as 
opposed  to  13.  5  sites  for  1971  data.  A  histogram  of  the  number  of  sites  for 
beamforming  from  1972  is  shown  in  Figure  II-2. 

V  Most  data  problems  seemed  to  be  associated  either  with 

"c^ad"  or  inverted  channels  (presumably  arising  from  renovation  of  seis¬ 
mometer  vaults)  or  tape  errors  incurred  during  recording  of  the  data  at 
SDAC  (mostly  in  late  July  and  early  August  1972).  Time  gaps  from  trans¬ 
mission  loss  or  NDPC  downtime  were  very  minor  problems. 

Processing  of  the  1971  events  was  done  using  three  different 
passbands  in  an  effort  to  determine  which  was  the  most  suitable.  It  was 
decided  that  the  passband  of  0.025-0.059  Hz  or  17-40  seconds  would  be  best 
and  was  used  for  processing  all  of  the  1972  events.  This  particular  band  was 
not  used  for  the  1971  events,  so  1971  results  from  the  0.020-0.059  Hz  band 
were  used  in  Table  II- 1. 
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FIGURE  II -1 

SITE  DIAGRAM  OF  THE  NORSAR  LONG-PERIOD  ARRAY 
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FIGURE  II- 2 

NUMBER  OF  SFTES  INCLUDED  IN  ARRAY  BEAM 
1972  EVENTS 


SECTION  III 


NOISE  ANALYSIS 


A.  INTRODUCTION 

The  results  presented  here  from  the  study  of  long-period 
noise  at  NORSAR  are  a  direct  continuation  of  the  work  reported  on  in 
Special  Report  No.  5.  By  combining  the  1971  data  reported  previously  with 
the  1972  data,  approximately  twenty  months  of  noise  history  have  been  analyzed. 
This  extensive  data  base  has  revealed  the  general  features  of  the  annual 
and  seasonal  variations  as  well  as  some  short-term  characteristics  of  the 
NORSAR  noise  field.  The  previous  study  also  suggested  some  areas  requiring 
a  closer  investigation,  such  as  the  correlation  of  noise  azimuths  to  meteor¬ 
ological  data  and  the  long  term  changes  in  noise  azimuth  and  RMS  level 
from  summer  to  winter.  These  phenomena  have  been  further  investigated 
and  detailed  results  are  presented  in  this  report. 

B.  DATA 

The  1972  data  used  in  this  study  are  listed  in  Table  HI-1. 

The  data  for  1971,  also  used  in  this  report,  are  listed  in  Special  Report  No.  5. 

The  1972  noise  samples  through  1  August  1972,  with  the 
exception  of  18  February  197  2,  wrre  several  hours  long.  The  samples  from 
18  August  1972  to  26  December  1972  were  42  minutes-40  seconds  long  and  were 
used  primarily  to  compute  average  RMS  amplitudes  to  the  end  of  1972.  The 
criteria  for  chosing  these  samples  are  given  in  Section  II. 

The  1972  data  were  processed  in  the  same  manner  as  the 
1971  data;  using  2.  0  second  sample  period  and  256  second  segments, 
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TABLE  III-l 


1972  NORSAR  NOISE  SAMPLES 


Day 

Start  Time 

hr:min 

Length 
hr: min:  sec 

Number  of 

available  sites 

5  January 

13:00 

6:24:00 

14 

5  January 

19:00 

6:24:00 

14 

1 3  January 

8:  30 

3:35:00 

13 

1  9  January 

4:  30 

4:45:52 

14 

30  January 

6:00 

5:49:52 

18 

1 0  F ebruary 

18:00 

4:  37:30 

19 

1 8  F ebruary 

5:  30 

0:21:08 

20 

26  F ebruary 

7:  30 

6:24:00 

15 

8  March 

22:00 

2:54:56 

16 

28  March 

22:30 

7:06:40 

18 

8  April 

10:  30 

5:02:56 

19 

23  April 

6:00 

4:45:52 

22 

1  7  May 

5:00 

1:  25:  20 

16 

27  May 

15:00 

3:07:44 

21 

3  June 

21:00 

4:1  6:00 

17 

1  5  June 

17:40 

3:35:00 

21 

7  July 

21:00 

2:42:08 

16 

1  3  July 

3:  20 

2:  50: 40 

20 

23  July 

14:  30 

2:29:  30 

21 

1  August  ’ 

14-  15 

3:0  3:28 

16 

1  8  August 

20:00 

0:42:40 

19 

28  August 

1:00 

0:42:40 

21 

6  September 

10:00 

0:42:40 

18 

1  7  September 

8:00 

0:42:40 

19 

28  September 

14:00 

0:42:40 

19 

1  7  October 

10:00 

0:42:40 

19 

27  October 

22:00 

0:42:40 

19 

6  November 

8:00 

0:42:40 

20 

1  5  November 

2  3:00 

0:42:40 

21 

26  November 

20:00 

0:42:40 

21 

6  December 

17:00 

0:42:40 

21 

26  December 

15:00 

0:42:40 

21 
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cross-power  matrices  were  formed  by  Fourier  transforming  segments, 
lianning,  crossmultiplying  the  transforms  and  stacking  the  segment  cross¬ 
power  matrices.  Dead  or  bad  sites  and  segments  were  eliminated  by  exam¬ 
ining  both  the  time  traces  and  the  power  spectra. 

The  data  were  not  corrected  for  instrument  response,  but 
the  nominal  amplitude  quantization  value  of  2.  47  m^/computer  counts  at 
0.  04  Hz  was  used  to  normalize  the  spectra. 

Only  the  vertical  components  were  used  in  this  analysis,  for 
the  following  reasons: 

•  The  1971  study  indicated  that  the  horizontals  generally  had 
similar  noise  levels  but  less  coherence  than  the  verticals 
in  the  20  to  40  second  band. 

•  The  requirement  that  all  three  components  be  operating  and 
free  of  non-seismic  noise  substantially  reduces  the  number 
of  sites  usable  for  f-k  spectral  analysis. 

•  A  visual  inspection  of  the  power. spectra  throughout  1972 
showed  that  the  horizontal -vertical  relationship  in  the  20 
to  40  second  band  had  not  changed  from  1971. 

Thus,  the  inclusion  of  the  horizontals  would  add  little  to  the  information  to 
be  derived  from  the  vertical  components. 

C.  SPATIAL  VARIABILITY 

Special  Report  No.  5  showed  that  the  NORSAR  array  was 
amplitude  equalized  within  +  2  dB  in  the  20  to  40  second  band  from  site  to  site 
during  1971.  Visual  inspection  of  the  broadband  power  levels  (0.00  to  0.  25  Hz) 

throughout  1972  and  a  check  of  the  1  August  1972  noise  sample  in  the  20  to  40 
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second  band  showed  that  amplitude  equalization  among  sites  continues  to 
be  within  ±  2  dB  during  1972. 


D.  TIME  VARIABILITY 

The  seasonal  dependency  of  the  noise  level  at  NORSAR  was 
determined  by  using  the  average  vertical  RMS  amplitude  in  the  20  to  40  second 
band.  This  average  was  computed  over  the  duration  of  each  noise  sample  for 
each  site  and  averaged  over  sites.  The  average  vertical  RMS  noise  level 
plotted  versus  day  of  the  year  is  shown  in  Figure  III — 1 ,  which  includes  1971 
data. 

The  primary  feature  of  Figure  III- 1  is  the  seasonal  variation 
in  noise  level.  The  noise  level  is  nearly  constant  from  1  May  1971,  when  data 
were  first  available,  to  27  November  1971  and  then  begins  to  rise  and,  while 
erratic,  stays  high  until  returning  to  a  nearly  constant  level  again  at  10  Februar 
1972.  Again  on  17  October  1972  the  RMS  level  begins  another  erratic  climb  to 
the  year's  end.  Thus  it  appears  that  the  period  between  February  or  March 
and  October  or  November  has  low  noise  levels  at  NORSAR,  with  average  RMS 
values  being  6.  8  m fj.  in  1971  and  7.  1  m fi  in  1972. 

An  investigation  of  the  short  term  noise  variations  in  the 
first  80  days  of  1972  was  made  from  the  noise  data  edited  with  the  signals 
processed  during  this  time.  The  noise  data  were  typically  only  1000  seconds 
(approximately  20  minutes)  long  and  preceded  the  P  arrival  on  site  1  vertical. 
The  data  were  bandpass  filtered  (20-40  second  band)  and  the  RMS  value  then 
computed.  Figure  III  -  2  shows  the  vertical  RMS  level  versus  day.  The  general 
structure  of  Figure  III- 2  agrees  with  that  of  Figure  III  - 1  in  that  maximum  and 
minimum  occur  on  the  same  day.  However,  the  general  'level  of  Figure  III- 2 
is  higher  than  Figure  HI-1  particularly  during  the  latter  part  of  February  and 
during  March.  Part  of  this  may  be  due  to  the  short  samples  of  Figure  IH-2 
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containing  some  low  level  coda  from  earlier  units  that  would  be  averaged  out 
over  the  longer  samples  of  Figure  HI-1.  Figure  III  -  2  illustrates  that  large 
short-term  variations  in  the  noise  level  occur  during  this  period;  clearly  the 
NORSAR  detection  capability  can  vary  by  up  to  about  0.  5  Mg  units  depending 
on  the  time  of  occurence  of  a  given  signal. 

The  very  high  noise  level  on  14  December  1971  was  investi¬ 
gated  in  more  detail  to  determine  the  duration  and  character  of  the  noise  im¬ 
mediately  preceding  and  following  this  noise  sample.  Data  were  edited,  when 
available,  from  1047Zof  13  December  (day  347)  to  0500Zof  15  December. 

The  broadband  (0  to  0.  25  Hz)  RMS  level  was  computed  over  successive  256 

second  gates.  Figure  III-3  shows  a  plot  of  the  RMS  level  of  every  third  gate 
versus  time. 

In  Figure  III  -  3 ,  noise  levels  are  fairly  constant  from  1100Z 
to  1 300 Z  on  13  December,  but  after  the  4  hour  data  gap  the  noise  begins  to  in¬ 
crease.  A  noise  peak  occurs  at  1100Z  on  14  December  and  then  gradually  the 
noise  level  drops  until  the  data  are  terminated  at  0400  Z0f  the  following  day. 

(A  swarm  of  Kamchatka  events  on  1  5  December  starting  at  0800Z  prevented 
editting  of  additional  noise  data.) 

Over  this  43  hour  period,  the  broadband  RMS  level  increased 
by  a  factor  of  10.  The  noise  in  the  20  to  40  second  band  showed  almost  as  much 
change,  with  RMS  amplitudes  of  16.  6  m/a  at  1200Z,  13  December;  77.  2m /a  at 
0500Z,  14  December  and  28.  2  m/a  at  0400Z,  15  December.  The  large  difference 
between  the  broadband  RMS  level  and  the  20-40  second  RMS  level  is  due  to  the 
large  6  second  peak  also  in  the  spectrum  during  this  time.  The  variation  of 
the  20-40  second  band,  above,  represents  a  change  of  0.  65  M  units  over  this 

g 

43  hour  period  and  represents  a  maximum  change  of  1.  2  M  units  relative  to 
the  NORSAR  minimum  of  6.  8  m/a  . 
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SITE  4  VERTICAL  BROADBAND  RMS  AMPLITUDE 
DURING  DECEMBER  13,  14,  AND  15,1971 


Figure  III- 3  also  shows  approximate  arrival  times  of  surface 
waves  from  five  events  listed  by  NOAA-PDE  with  mb>  4.  5  and  occuring  in 
this  time  period.  These  events  were  not  routinely  processed  so  that  beam  de¬ 
tection  data  were  not  available,  however,  on  the  raw  time  traces  only  the  Masca 
rene  Is.  event  arriving  at  1047Z,  14  December  was  definitely  visible.  The 
others  are  not  visible  probably  due  to  their  depths. 

The  data  over  this  43  hour  period  supplement  the  data  in 
Figure  III-2  in  that  while  the  day-to-day  changes  in  RMS  level  may  be  large, 
the  hourly  change  is  gradual  and  not  as  sudden  as  Figure  III- 2  might  suggest. 

E.  NOISE  DIRECTIONALITY 

Special  Report  No.  5  showed  that  the  noise  direction  on  14 
December  appeared  to  be  directly  correlated  with  areas  of  low  pressure  and 
large  swell  heights  in  the  North  Atlantic.  From  the  data  used  in  Figure  III-3, 
high  resolution  f-k  spectra  using  maximum  likelihood  filters  (Barnard,  1969) 
were  generated  at  0.  06  Hz  for  the  time  intervals  from  1047Z  to  1350Z,  13 
December;  0400Z  to  0603Z,  14  December;  and  0440Z  to  0513Z,  14  December. 
These  spectra,  shown  in  Figure  III-4,  indicate  a  noise  peak  at  3.  5  km/sec. 
which  shifts  clockwise  from  280°  to  350°  to  10°  azimuth  over  this  43  hour  period. 

During  this  time  period,  National  Weather  Service  surface 
pressure  charts  and  Combined  Sea  and  Swell  Height  Charts  were  available. 
Figure  III- 5  is  a  map  of  the  North  Atlantic  <"cean  region  on  which  the  move¬ 
ments  of  the  three  low  pressure  areas  present  on  these  dates  are  plotted. 

The  two  areas  of  maximum  wave  heights  (taken  from  the  two  wave  height 
charts  available)  also  are  shown. 

During  this  43  hour  period,  the  northernmost  low  pressure 
area  (labeled  C-C')  and  the  northernmost  area  of  maximum  wave  height  move 
from  west  to  north  relative  to  NORSAR  coincident  with  the  west  to  north  move¬ 
ment  of  the  micros eismic  peak  azimuth  at  NORSAR.  This  indicates  that  the 
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ber  13  at  1000Z 


FIGURE  III- 5 

PATHS  OF  LOW  PRESSURE  CENTERS  IN  6  HOUR  INTERVALS 
AND  MAXIMUM  WAVE  HEIGHTS  13-15  DECEMBER  1971 
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observed  noise  azimuth  is  directly  correlated  with  either  the  low  pressure 
area  or  the  area  of  maximum  wave  height. 

The  suggestion  of  a  direct  correlation  between  noise  azimuth 
and  either  maximum  wave  height  area  or  low  pressure  area  was  further  in¬ 
vestigated  in  the  noise  of  5  January  1972.  The  noise  sample  of  5  January 
consists  of  1 1  hours  of  noise  starting  at  1  300Z  of  5  January  (omitting  the  data 
between  1  500  Z  -  1600Z  which  contained  an  event).  High  resolution  f-k  spectra 
at  the  microseismic  peak  frequency  were  generated  hourly  and  the  results  by 
hour,  frequency  and  peak  azimuth  are  tabulated  in  Table  III  -  2.  The  velocities 
of  all  peaks  were  near  3.  5  km/ sec  .  The  significant  feature  of  Table  III-2 
is  the  absence  of  a  continuous  range  of  azimuths  but  rather  a  grouping  of  a- 
zimuths  from  285°-  290°  and  3  50°  to  0°  (the  same  as  the  13-15  December 
period).  One  value  at  302  is  an  exception. 

Low  pressure  centers  and  maximum  wave  height  areas  from 
the  National  Weather  Service  charts  for  this  period  are  plotted  in  Figure  III-6. 
The  two  major  azimuths  of  285°  and  0°  correlate  very  well  with  the  low  pres¬ 
sure  centers  although  the  northern  low  (labeled  B-B')  is  very  weak.  The 
single  peak  at  302°  could  be  due  to  the  high  seas  off  Iceland. 

Of  additional  interest  is  the  structure  of  the  individual  high 
resolution  f-k  spectra  of  these  samples.  Figure  III  -  7  is  a  series  of  plots  of 
the  power  values  in  the  high  resolution  f-k  spectra  at  a  velocity  of  3.  5  km/ sec 
for  several  of  the  5  January  samples.  The  power  values  are  in  dB  relative  to 
the  peak  value  and  plotted  versus  azimuth.  These  data  show  ihree  local  max¬ 
imum  in  the  spectra,  at  200°to  220°,  280°  tok  300°,  and  0°  azimuth.  Figure 
III  -  7  suggests  that  the  above  three  azimuthal  noise  directions  are  always  pre¬ 
sent  in  the  noise  spectral  but  at  varying  power  levels. 


TABLE  III- 2 


HOURLY  PEAK  FREQUENCY  AND  PRIMARY  SOURCE  AZIMUTH 
OF  THE  MICROSEISMIC  NOISE  ON  5  JANUARY  1972 


Hour  (GMT) 
(Approximate) 

F  requency 
(Hz) 

Azimuth 

(Deg) 

1300-1400 

0.  066 

285° 

1400-1500 

0.066 

o 

285 

1600-1700 

0.  066 

o 

290 

1700-1800 

0.  066 

o 

0 

1800-1900 

0.  066 

355° 

1900-2000 

0.062 

350° 

2000-2100 

0.  062 

30  2° 

2100-2200 

0.  066 

o 

358 

2200-2300 

0. 066 

285° 

2300-0030 

0.066 

285° 
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FIGURE  III  -  6 

PATHS  OF  LOW  PRESSURE  CENTERS  IN  SIX  HOUR  INTERVALS 
AND  MAXIMUM  WAVE  HEIGHTS  5  JANUARY  1972 
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Azimuth  (Deg) 


FIGURE  III -7 


AZIMUTHAL  POWER  DISTRIBUTION  AT  3.  5  KM/SEC. 
FROM  1400  Z  TO  1900  Z  ON  5  JANUARY  1972 
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The  other  noise  samples  of  Table  III- 1  also  *end  to  reinforce 


the  suggestion  of  discrete  noise  azimuths  rather  than  a  continuous  range. 

Table  III  -  3  is  a  list  of  the  1972  long  noise  samples  with  the  high-resolution 

f-k  spectra  power  peaks  listed  by  frequency  of  the  peak  and  azimuth  (all 
the  velocities  were  at  3.5  km/sec).  From  the  data  of  Table  III  -  3  and  the 
1971  data  the  histogram  of  Figure  III  -  8  was  constructed  with  the  number  of 
occurences  of  peak  power  tabulated  in  10°  azimuthal  increments.  This  histo¬ 
gram  is  further  divided  into  shaded  and  unshaded  samples.  The  shaded  samples 
are  data  from  15  May  through  8  October  1971  and  13  July  through  1  August  1972 
and  contain  all  the  easterly  noise  peak  azimuths.  It  is  of  interest  to  note  that 
the  above  time  gates  for  easterly  noise  are  not  the  same  as  the  times  of  mini¬ 
mum  noise  at  NORSAR  so  that  the  lowering  of  the  noise  level  does  not  coincide 
with  the  beginning  of  easterly  peak  noise  azimuths.  The  above  time  gates  for 
easterly  noise  azimuths  nearly  correspond  with  the  summer  season  and  are 
termed  the  summer  noise;  data  outside  these  times  are  called  winter  noise. 

In  the  histogram  of  Figure  III  —  8  there  are  five  azimuths 

which  contain  approximately  70%  of  the  noise  peaks;  0°,  60°,  90°,  250°  and 

o  n 

280  .  Of  the  winter  samples  80%  of  the  noise  peaks  are  at  azimuths  of  0  , 

250  and  280°. 

Figure  III-9  shows  relative  power  versus  azimuth  taken  from 
the  high-resolution  f-k  spectra  at  3.  5  km/ sec  for  four  winter  samples.  In 
all  the  examples  shown,  and  in  all  the  winter  noise  f-k  spectra  investigated 
one  or  more  of  the  previously  noted  azimuths  appear  as  primary  or  secondary 
maxima.  The  fact  that  these  azimuths  consistently  appear  in  a  wide  variety 
of  noise  samples  and  winter  weather  conditions  gives  strong  evidence  that  the 
winter  noise  has  only  a  few  discrete  azimuthal  directions  0°,  250°  and  280°. 

•  The  existence  of  these  discrete  azimuths  are  of  significant 
importance  in  the  generation  of  long  period  noise.  Because  the  low  pressure 
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TABLE  III -3 


VERTICAL  PEAK  POWER  AZIMUTH, FREQUENCY, 
AND  VELOCITY  BY  DAY  OF  YEAR  1972 


Date  (1972) 

Frequency 

(Hz) 

Peak  Azimuth 
(Deg.) 

l  2  January 

.  066 

0° 

1  9  J  anuary 

.  063 

300° 

JO  January 

.  063 

283° 

10  F ebruary 

.  0  66 

255° 

26  F ebruary 

.066 

o 

250 

8  March 

.066 

o 

5 

28  March 

.  070 

o 

285 

8  April 

.  078 

218° 

23  April 

.063 

285° 

17  May 

.078 

0° 

27  May 

.0  70 

255° 

3  June 

.074 

265° 

1 5  June 

.055 

335° 

7  July 

.  070 

255° 

13  July 

.  055 

93° 

23  July 

.  066 

65° 

1  August 

.059 

95° 

Azimuth  (Deg) 


FIGURE  III -8 

HISTOGRAM  OF  VERTICAL  PEAK  POWER  AZIMUTHS 
FROM  NOISE  SAMPLES  IN  1971  AND  1972 


areas  and  maximum  wave  height  areas  move  continuously  northeastward  in 
the  North  Atlantic  a  concentration  of  energy  coupling  at  various  points  along 
the  Norwegian  coast  due  to  ocean  bathymetry  appears  to  be  the  probable 
factor.  This  effect  is  to  be  further  investigated  in  future  wo.-k. 

As  stated  previously  the  summer  noise  samples  contain 
all  the  easterly  noise  peak  azimuths  (Figure  IH-8).  Even  though  the  summer 
peaks  show  a  predominantly  easterly  direction  the  noise  is  generally  almost 
isotropic.  Figure  HI-10  is  a  contoured  plot  of  the  high-resolution  f-k  spectra 

of  23  July  1972.  The  noise  in  Figure  III-10  is  nearly  isotropic  between  3.  5 
and  4.  0  km/ sec. 

Note  that  the  distinction  between  winter  and  summer  noise  was 
made  on  the  basis  of  noise  azimuth  rather  than  noise  level.  The  lowering  of 
the  RMS  noise  level  to  the  base  levels  shown  in  Figure  IH-1  does  not  coincide 
with  the  change  to  easterly  azimuths.  In  Figure  III-l  the  RMS  level  at  10 
February  1972  is  near  its  minimum  value  of  7.  1  m fx  but  the  structure  of  the 
high -resolution  f-k  spectra  on  this  day  (Figure  IH-9)  shows  the  strong  dis¬ 
crete  noise  azimuths  of  winter.  The  low  RMS  noise  level  on  10  February 
with  the  absence  of  an  easterly  noise  peak  suggests  that  the  summer  noise,  at 
a  comparable  RMS  to  10  February,  can  not  be  a  constant  year  around  back¬ 
ground  level  but  is  due  to  a  different  generating  mechanism  than  that  of  the 

winter  noise,  .lowever,  the  source  mechanism  for  summer  noise  is  presently 
unknown. 

In  summary,  the  following.conclusions  from  this  noise  analysis 

can  be  made: 

•  Site-to-site  variation  continues  to  be  within  +  2dB  in  the  20-40 

second  band  for  1972. 
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RMS  noise  levels  are  seasonally  dependent.  The  first  and 
last  two  months  of  each  year  are  characterized  by  high  noise 
levels  and  large  changes,  often  occuring  daily,  whereas 
the  remainder  of  the  year  has  noise  levels  uniformly  low. 

Noise  levels  appear  as  a  function  of  the  severity  of  North 
Atlantic  weather  conditions. 

Winter  noise  generally  appears  to  be  from  one  or  more  of 
three  azimuths  0  ,  250  and  280  .  The  reason  for  these 
preferred  azimuths  will  be  further  investigated. 

Summer  noise  is  nearly  isotropic  with  a  weak  maximum 
to  the  east. 

Seasonal  changes  in  noise  azimuth  are  not  i  eflected  in 
concurrent  changes  in  noise  level.  Noise  azimuths  show 
little  correlation  with  RMS  level. 

Winter  detection  capabilities  can  change  by  0.  50  -  0.  65M 

s 

units  from  day  to  day  and  can  be  up  to  1. 2M  units  higher 
than  those  observed  during  the  summer. 


Ill  -  22 


SECTION  IV 

ARRAY  PROCESSING  PERFORMANCE 


A.  INTRODUCTION 

The  effectiveness  of  array  processing  at  NORSAR  was 
evaluated  by  comparing  the  performances  of  a  conventional  beamsteer  (BS) 
processor  and  a  multichannel  Wiener  filter  (MCF)  processor.  Results  from  last 

year's  work,  discussed  in  Special  Report  No.  5.  showed  that,  under  winter 

noise  conditions.  MCF  processing  gave  significantly  better  noise  reduction 
than  beamsteering.  The  purpose  of  the  evaluation  this  year  was  to  see  if 
MCF's  continued  to  significantly  outperform  beamsteering  and  to  explore 
the  feasibility  of  routinely  designing  and  applying  multichannel  filters.  The 
following  analyses  were  undertaken  to  accomplish  these  objectives: 

Estimation  of  the  noise  rejection  achieved  by  MCF  and  BS 
processors  for  various  array  and  subarray  configurations. 

Examination  of  the  azimuthal  dependence  of  array  gains. 

Comparison  of  signal -to -noise  ratio  (SNR)  improvements 

obtained  by  MCF  and  BS  processors  for  several  seismic 
events . 

Multichannel  filters  and  beamsteer  processors  were  designed 
from  and  applied  to  eleven  winter  noise  samples  to  measure  noise  reduction 
for  four  array  and  subarray  configurations:  the  full  array,  the  outer  ring 
sites,  the  inner  ring  sites,  and  a  special  subarray  of  seven  close  sites.  One 
of  the  noise  samples  was  used  to  study  the  azimuthal  dependence  of  array  gains. 


IV-1 


In  order  to  test  the  effects  of  these  processors  on  signal  waveforms,  MCF 
and  BS  processors  were  also  designed  and  applied  to  four  events  (and  their 
associated  noise  samples). 

B.  WINTER  NOISE  REJECTION  BY  MCF  AND  BS  PROCESSORS 

The  present  work  on  array  processing  performance  was 
directed  toward  the  noise  rejection  achieved  by  MCI  and  BS  processors  on 
winter  noise  at  NORSAR,  using  various  array  and  subarray  configurations. 
Eleven  typical  winter  samples,  listed  in  Table  IV-1,  were  used  for  the  noise 
rejection  analysis.  These  samples  covered  the  period  from  5  January  1972  to 
8  March  1972.  Note  that  for  the  last  two  of  these  samples,  noise  levels  had 
returned  to  "summer"  levels,  although  noise  azimuths  were  still  typical  of 
wintertime  conditions  (see  Section  III).  The  shortest  noise  sample  was  five 
hours  long  (71  data  segments)  and  the  longest  extended  approximately  7  hours 
(100  segments).  The  time  data  were  checked  for  quality  and  for  surface  waves 
from  unreported  events.  The  data  from  bad  sites  and  the  portions  of  the  data 
including  signals  or  glitches  were  deleted. 

The  MCF's  were  designed  for  only  the  vertical  components 
of  the  Rayleigh  wave.  The  following  MCF  design  parameters  were  used  '.o 
allow  comparison  of  results  to  previous  NORSAR  and  ALP  A  results. 

•  Frequency-domain  design  by  transforming,  cross-multiplying, 
and  stacking  256-second  segments  at  all  frequencies  (0.0  - 
0.25  Hz),  The  cross- spectra  were  not  smoothed. 

•  Dispersive  signal  model  oriented  to  a  beam  direction  of  90°, 
the  approximate  azimuth  to  central  Asia, 

•  Signal-to-noise  ratio  equal  to  four  at  all  frequencies,  to  main¬ 
tain  a  balance  between  signal  preservation  and  noise  rejection. 

•  Two  percent  white  noise  added  to  the  data  to  stabilize  filter 
design. 
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NOISE  SAMPLE  DATA  USED  TO  ESTIMATE  ARRAY  PROCESSING  PERFORMANCE 


£  1  segment  is  256  seconds  long 


The  frequency  domain  filter  coefficients  were  transformed  to  time  domain 
operators  and  subsequently  were  convolved  with  time  data  to  obtain  a  MCF 
output.  The  MCF  and  BS  processors  were  applied  to  two  sections  (each 
about  3000  seconds  long)  of  each  noise  sample.  The  first  section,  within 
the  MCF  design  gate,  was  called  the  "on-design-  portion;  the  second  section, 
outside  of  the  MCF  design  gate  but  within  the  same  noise  sample,  produced 
the  "off-design-  protion.  Table  IV-2  gives  a  list  of  the  sites  used  (see  also 
Figure  II-l  in  Section  II  of  this  report). 

Bandpass  filters  then  were  applied  to  the  beams  and  the 
RMS  values  were  calculated.  The  passbands  were  0.025  -  0.050  Hz,  0.025 
0. 059  Hz,  and  0.  020  -  0.  10  Hz.  The  first  passband  has  been  used  for  general 
noise  analysis  and  was  used  here  for  comparison;  the  second  passband  is 
normally  used  for  event  processing  at  NORSAR;  the  third  passband  covers 
the  central  portion  of  the  NORSAR  long-period  seismometer  response. 

The  noise  reduction,  i.  e.,  array  gain,  achieved  by  each  pro¬ 
cessor  was  computed  as  the  ratio  of  average  single  site  RMS  amplitude  to 
the  RMS  amplitude  of  the  beam.  The  MCF  over  BS  noise  reduction  is  listed 
in  Table  IV-3  for  various  array  configurations.  Within  the  on-design  portion 
of  the  data,  the  MCF  was  alwnys  superior  to  the  beams  teer.  This  superior¬ 
ity  of  the  MCF  was  generally  maintained  with  data  in  the  off-design  gate, 
although  improvements  were  smaller.  The  broadest  band  (0.  02  -  0. 10  Hz) 
MCF  had  the  best  noise  reduction  and  the  narrowest  band  (0.  025  -  0.  050  Hz) 
MCF  showed  the  poorest.  This  is  expected  because  the  coherent  energy  is 
concentrated  between  0.05  and  0.10  Hz. 

The  correlation  between  the  noise  level  and  the  MCF  over 
BS  noise  reduction  is  worth  noting.  Using  the  full  array  results  (Table  IV-3), 
the  average  single  site  RMS  noise  values,  the  beamsteer  levels  for  off-design 


SITES  USED  FOR  ARRAY  PROCESSING  EVALUATION 
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LIST  OF  SITES  USED  FOR  ARRAY  PROCESSING  EVALUATION 
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TABLE  IV -3 


IMPROVEMENTS  (dB)  IN  NOISE  REDUCTION 
OF  MCF  OVER  BS  PROCESSOR 
(PAGE  1  OF  2) 


Noise 

Sample 

Site 

Configuration 

Number 

of 

Sites 

MCF  over  BS  Noise  Reduction 
(in  dB) 

On-Design 

Off-Design 

0.  025 

to 

0.  050 

Hz 

0.  025 

to 

0.  059 

Hz 

0.  020 

to 

0.  10 

Hz 

0.  025 

to 

0.  050 

Hz 

0.  025 

to 

0.  059 

Hz 

0.  020 

to 

0.  10 

Hz 

NOIS 

Full  Array 

12 

1.  8 

3.  5 

7.  2 

1.  0 

0.7 

1.  2 

*  0  0  5 

Outer  Ring 

9 

2.  9 

2.  9 

5.  9 

2.  8 

4.2 

6.  1 

*13  A 

Inner  Ring 

5 

2.2 

3.  5 

5.  5 

’  5 

(-0.3) 

0.  2 

Special  Subarray 

7 

2.3 

2.  7 

7.  5 

Z .  4 

2.  1 

(-1.2) 

NOIS 

Full  Array 

11 

5.  1 

3.  5 

6.6 

3.  1 

2.  3 

8.  5 

*005 

Outer  Ring 

8 

3.  5 

6.  0 

4.  4 

2.  9 

3.  9 

4.6 

*19  A 

Inner  Ring 

6 

7.  5 

6.  5 

6.  1 

7.  4 

5.  6 

4.  3 

Special  Subarray 

7 

4.  0 

4.  2 

6.  4 

6.  0 

6.4 

6.  3 

NOIS 

Full  Array 

14 

3.  0 

4.  3 

7.  6 

0.  2 

0.  7 

2.  7 

*012 

Outer  Ring 

8 

4.1 

4.  7 

6.  3 

1.  2 

2.  0 

4.  1 

*21 A 

Inner  Ring 

7 

3.  8 

5.  0 

5.  8 

0.  6 

1.2 

2.  3 

Special  Subarray 

7 

4.  5 

5.  0 

5.  8 

0.  6 

0.  5 

2.1 

Full  Array 

15 

2.  6 

3.  0 

4.  0 

0.  0 

0.  5 

3.  3 

■lllii 

Outer  Ring 

9 

2  3 

2.  9 

3.  7 

0.  2 

1.4 

3.  7 

*03A 

Inner  Ring 

7 

2.  1 

2.  3 

3. 1 

0.  6 

0.  7 

2.  7 

Special  Subarray 

6 

3.  8 

3.  3 

4.  5 

(-0.1) 

(-0.1) 

5.  0 

NOIS  | 

Full  Array 

13 

8.  1 

7.  8 

8.  4 

6.  2 

5.  9 

6.  5 

*013 

Outer  Ring 

8 

3.  5 

4.  5 

5.  0 

3. 1 

3.  3 

4.  0 

*08A 

Inner  Ring 

5 

2.  1 

3.  8 

4.  3 

2.4 

3.  9 

4.  0 

Special  Subarray 

7 

2.  6 

3.  1 

5.  4 

1 .  3 

2.  5 

5.  3 
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IMPROVEMENTS  (dB)  IN  NOISE  REDUCTION 
OF  MCF  OVER  BS  PROCESSOR 
(PAGE  2  OF  2) 


Noise 

Sample 

Site 

Configuration 

Number 

of 

Sites 

MCF  over  BS  Noise  Reduction 
(in  dB) 

On-Design 

Off-Design  i 

0.  025 

to 

0.  050 

Hz 

0.  025 

to 

0.  059 

Hz 

0.  020 

to 

0.  10 

Hz 

0.  025 

to 

0.  050 

Hz 

0.  025 

to 

0.  059 
Hz 

0.  020 

to 

0. 10 

Hz 

NOIS 

Full  Array 

13 

1.  6 

5.  0 

(-1.1) 

B 

2.  9 

*019 

Outer  Ring 

9 

2.  1 

4.  3 

0.  3 

m 

6.4 

*04  A 

Inner  Ring 

6 

1.  5 

3.4 

1.  2 

1.2 

3.  1 

Special  Subarray 

7 

2.  5 

2.  8 

0.  8 

0.  8 

3.  4 

NOIS 

F ull  Array 

12 

3.  1 

6.  6 

3.2 

HI 

*019 

Outer  Ring 

8 

3.6 

6.  0 

6.  5 

1.  3 

4.1 

*10  A 

Inner  Ring 

5 

1.  8 

3.2 

6.3 

1.  7 

1.  6 

i 

Special  Subarray 

7 

1.  2 

2.  5 

6.6 

0.  6 

1.  0 

NOIS 

Full  Array 

15 

2.  5 

6. 1 

. 

10.  3 

2.  4 

3.1 

*030 

Outer  Ring 

10 

2.  5 

5.  2 

8.  1 

2.  2 

3.6 

5.  8 

*2  3  A 

Inner  Ring 

6 

3.  9 

7.  0 

9.  8 

5.4 

7.  8 

Special  Subarray 

7 

3.1 

5.6 

8.6 

2.  6 

4.  0 

8.2 

NOIS 

Full  Array 

14 

3.  6 

m 

8.  5 

2.4 

4.  9 

9.  0 

*041 

Outer  Ring 

8 

1.  5 

4.2 

0.  9 

2.2 

3.  5 

*18  A 

Inner  Ring 

6 

5.  3 

6.1 

S.  0 

6.6 

5.  5 

Special  Subarray 

7 

5.  4 

kh 

7.  0 

4.6 

5.  8 

6.1 

NOIS 

F ull  Array 

14 

3.  6 

3.  8 

5.  6 

(-1.4) 

SO 

2.  9 

*057 

Outer  Ring 

8 

4.  6 

4.  0 

5.  2 

(-1.3) 

m 

3.  9 

*07  A 

Inner  Ring 

7 

5.  9 

4.  2 

4.  4 

2.  7 

2.  5 

3.  6 

Special  Subarray 

7 

6.  3 

5.  9 

6.  2 

2.  2 

3.  3 

NOIS 

F ull  Array 

16 

1.  9 

3. 1 

6.4 

0.4 

1.  8 

5.  6 

*068 

Outer  Ring 

3.  5 

4.  0 

4.  0. 

0.  7 

2.  8 

2.1 

*2  2  A 

Inner  Ring 

6 

3.  8 

3.4 

5.  5 

3.  9 

3.  8 

4.  9 

Special  Subarray 

7 

3.  8 

3.  9 

6.  5 

2.1 

. 

2.  2 

5.  9 

and  on-design  noise,  and  the  MCF  over  BS  noise  reduction  all  were  plotted 
versus  the  noise  samples  for  each  of  the  three  passbands  (Figures  IV-1  to 
IV-3).  For  off-design  noise,  the  MCF  improvements  range  from  -1  to  6  dB 
in  the  0.  025  -  0.  050  Hz  passband,  0  to  6  dB  in  the  0.025  -  0.059  Hz  passband, 
and  1  to  8  dB  in  the  0.  02  -  0. 10  Hz  passband.  The  on-design  improvements 
were  generally  better  than  the  off-design  improvements,  as  expected.  Sur¬ 
prisingly,  the  best  MCF  improvements  were  not  obtained  when  the  noise 
levels  were  the  highest. 

Table  IV-4  shows  the  mean  improvements  of  the  MCF  over 
BS  processor  for  several  array  and  subarray  configurations.  For  any  given 
passband  the  on-design  noise  always  had  better  improvements  than  the  off- 
design  noise.  For  both  off-design  and  on-design  portions,  the  results  indi¬ 
cate  that  the  inner- ring  and  the  seven- site  subarray  yielded  slightly  better 
MCF  over  BS  improvements  than  the  outer-ring  and  full  array.  Also,  the 
Improvements  listed  are  about  the  same  as  those  observed  for  last  year's 
data. 

C.  AZIMUTHAL  DEPENDENCE  OF  ARRAY  GAIN 

The  MCF  and  BS  processor  performances  would  be  expected 
to  be  strongly  dependent  on  noise  field  structure  relative  to  the  direction  of 
the  beam.  To  examine  this  aspect  the  noise  sample  NOIS*019*10A  was  ana¬ 
lyzed  in  detail.  The  high  resolution  f-k  spectrum  for  this  noise  sample 
showed  that  at  0.  050  Hz  and  0.  059  Hz  the  noise  came  from  northwest,  with 
a  peak  at  292°. 

The  MCF  and  BS  outputs  were  formed  using  12  sites  for  look 
directions  of  0°,  45°,  90°,  135°,  180°,  225°,  270°  and  315°.  These  beams 
included  both  on-design  and  off-design  portions  of  the  noise  sample.  The 
resulting  MCF  and  BS  array  gains  are  shown  in  Figures  IV-4  through  IV-9. 
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FIGURE  IV-1 

PROFILES  OF  ARRAY  GAINS 
(0.  025  -0.  050  Hz) 
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Beam  RMS  (m fi  )  MCF  over  BS  Improvements  (dB) 
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IMPROVEMENT  vs.  AZIMUTH  FOR  MCF  AND  BS  PROCESSORS 
0.  025-0.  050  Hz,  On-Desien  Noise,  Day  019,  1972 
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IMPROVEMENT  vs.  AZIMUTH  FOR  MCF  AND  BS  PROCESSORS 
0.  025-0.  059  Hz,  Off-Design  Noise,  Day  019,  1972 
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IMPROVEMENT  vs.  AZIMUTH  FOR  MCF  AND  BS  PROCESSORS 
0.  020-0.  10  Hz,  Off-Design  Noise,  Day  019,  1972 


From  Figure  IV-4  to  Figure  IV-6  (on-design  noise  results),  the  MCF  had 

from  2  to  5  dB  better  noise  reduction  for  azimuth  from  0°  to  270°  and  essen- 

o  o 

tially  equivalent  performance  between  270  and  0  .  For  off-design  gains 
(Figures  IV-7  to  IV-9),  the  MCF's  still  did  better  than  BS  processors,  but 
the  gains  were  0  to  4  dB.  Improvements  for  both  MCF  and  beamsteer  pro¬ 
cessors  were  smaller  between  270°  and  360°,  which  would  be  expected  since 
that  is  the  sector  where  most  of  the  propagating  energy  is  coming  from. 

D.  MCF  AND  BS  PROCESSORS  APPLIED  TO  SIGNALS 

Selected  array  processors  were  also  applied  to  a  few  signals 
to  compare  their  performances  for  practical  signal  extraction  problems. 

Care  ^vas  taken  to  insure  that  channels  used  in  the  design  of  the  MCF  also 
were  available  for  event  processing.  Time  separations  between  signals  and 
noise  samples  ranged  from  8  hours  to  31  hours. 

Noise  rejection  estimates  were  obtained  by  computing  RMS 
levels  in  a  1500-second  time  gate  taken  prior  to  P-wave  arrival.  Relative 
signal  beamforming  loss  was  estimated  by  measuring  the  maximum  zero- 
to-peak  values  of  the  signal  waveforms.  Results  are  presented  in  Table  IV-5. 
The  MCF  signal  amplitudes  were  essentially  the  same  as  those  of  the  beam- 
steer,  indicating  that  signal  loss  from  beamforming  was  the  same  for  MCF  and 
BS  processors.  However,  somewhat  greater  noise  reductions  were  achieved 
by  the  MCF.  The  net  improvements  in  SNR  achieved  by  the  MCF  over  the 
BS  ranged  from  0  to  6  dB,  with  an  average  at  about  3  dB.  This  is  on  the 
same  order  as  off-design  mean  improvements  shown  in  Table  IV-4,  despite 
the  fact  that  these  noise  data  generally  were  further  separated  in  time 
from  the  design  noise  than  the  samples  used  in  Table  IV-4.  It  ap¬ 
pears  that  MCF  processing  can,  on  the  average,  produce  some  improvement 
over  BS  processing  for  winter  noise  provided  that  a  long  noise  sample  is 
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used  in  the  MCF  design  and  that  this  noise  sample  is  within  about  one  day 
of  the  event  to  be  processed. 

E.  SUMMARY 

The  main  results  from  the  comparative  evaluation  of  the 
MCF  and  BS  processor  using  winter  noise  are  the  following: 

•  MCF' s  generally  achieved  better  noise  rejection  than  beam- 
steer  filters:  2-8  dB  for  on-design  noise  and  0~6  dB  for 
off-design  noise  in  the  0.025  -  0.059  Hz  signal  processing 
band. 

•  The  MCF  outperformed  the  beamsteer  processor  by  2~5  dB 
for  on-design  noise  and  by  0~4  dB  for  off-design  noise  fov 
look  directions  significantly  different  from  the  predomin¬ 
ant  noise  azimuths.  Serious  deterioration  in  array  gains 
for  both  the  MCF  and  beamsteer  was  seen  at  azimuths 
spanning  predominant  noise  azimuths. 

•  The  MCF  gave  0~6  dB,  and  typically  3  dB  more  SNR  im¬ 
provement  than  the  BS  for  example  of  "practical"  signal 
extraction  (usually  real  signals),  suggesting  that  with  the 
aid  of  MCF  processing  the  winter  detection  threshold  at 
NORSAR  could  be  decreased  by  on  the  order  of  0.  15  M 
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MATCHED  FILTERING  PERFORMANCE 

A.  INTRODUCTION 

Matched  filters  were  applied  to  the  beams  formed  for  the  events 
processed  this  year  to  measure  the  improvement  in  signal-to-noise  ratio  (SNR) 
which  could  be  obtained  at  NORSAR.  The  objectives  of  measuring  the  effective¬ 
ness  of  matched  filters  were: 

•  To  determine  the  available  SNR  gain  from  matched  filters. 

•  To  examine  the  variation  in  these  gains  by  region  and  surface 
wave  mode. 

•  To  determine  if  seismic  regions  could  be  catagorized  by  matched 
filter  parameters. 

These  objectives  were  reached  with  good  success.  In  particular,  we  were 
highly  successful  in  our  first  attempt  to  define  seismic  regions  on  the  basis  of 
a  matched  filter  parameter. 

Chirp  filters  were  routinely  applied  to  the  three-component 
beams  formed  for  each  event  processed  this  year.  In  contrast  to  last  year's 
work,  reference  waveform  filter s* (RWF)  were  not  applied  routinely  to  all 
events  from  the  areas  of  central  Italy,  eastern  Kazakh,  and  a  small  area  east 
of  Kamchatka.  The  reasons  for  this  were  that: 


*  In  the  report  on  last  year's  work,  reference  waveform  filters 

were  called  "master"  waveform  filters. 
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•  RWF's  generally  did  not  show  a  clear-cut  superiority  to  chirp 
filters  in  signal-to-nois e  (SNR)  improvement. 

•  The  event  source  parameters,  which  control  the  effectiveness  of 
a  RWF,  were  unknown;  thue  there  was  no  a  priori  means  of 
selecting  a  RWF  which  would  produce  good  SNR  improvement. 

•  The  relative  importance  of  the  source  parameters  and  their 
degree  of  interaction  also  were  unknown;  hence  they  could  not 
be  used  to  categorize  event  behavior  on  a  regional  basis. 

The  size  of  the  combined  data  base  this  year  allowed  a  suffici¬ 
ently  dense  sampling  of  events  in  some  areas  so  that  distinct  patterns  of 
behavior  from  region  to  region  could  be  ascertained.  These  patterns  were 
primarily  due  to  similarities  in  chirp  filter  length  and  SNR  improvements  and 
permitted  tentative  identification  of  regional  boundaries.  The  events  in  each 
region  were  analyzed  in  terms  of  SNR  improvement. 

B.  DESCRIPTION  OF  METHOD 

The  SNR  improvement  of  the  matched  filter  was  computed  by 
applying  in  turn  the  matched  filter  and  a  bandpass  filter  to  the  signal  waveform 
and  to  a  section  of  noise  preceding  the  signal.  The  SNR  of  the  output  of  each 
filter  was  then  computed  on  the  ratio  of  zero-to-peak  signal  amplitude  to  RMS 
noise.  The  ratio,  in  dB,  of  the  matched  filter  SNR  to  the  bandpass  filter  SNR 
is  the  SNR  improvement.  The  signal  gate  generally  included  all  major  signal 
energy;  the  noise  gate  was  generally  1000  seconds  long  and  was  chosen  to 
exclude  any  obvious  signal  energy  from  other  events. 

% 

The  number  describing  peak  signal  amplitude  more  correctly 
describes  peak  signal-plus-noise  amplitude,  particularly  for  weak  signals. 

For  this  reason,  we  shall  henceforth  refer  to  the  signal-plus-noise-to-noise 
ratio  (SNNR)  instead  of  signal-to-noise  ratio. 


The  filters  were  designed  in  the  frequency  domain  by  sper'fying 
the  passband  of  the  bandpass  filter  and  the  passband  and  lengthf'of  the  chirp. 

The  responses  were  not  tapered  at  the  band  edges.  The  passband  of  the 
bandpass,  RWF ,  and  chirp  filters  i  sed  for  signal  processing  at  NORSAR  this 
year  was  0.  025-0.  059  Hz,  or  17-40  seconds.  The  chirps  were  linear  in  the 
sense  that  the  phase  rate  of  change,  i.  e.,  frequency,  was  a  linear  mnction  of 
time. 

The  beams  to  which  these  filters  were  applied  were  formed  using 
a  phase  velocity  of  3.  5  km/sec  for  the  vertical  and  radical  components  and 
4.0  km/sec  for  the  transverse  component. 

Five  chirp  filters  and  one  bandpass  filter  were  applied  sequentially 
to  each  component  beam.  Trial  chirp  lengths  were  chosen  to  c  over  a  moderate 
range  about  the  predicted  best  length,  which,  was  obtained  either  from  group 
velocity  curves  or  from  previous  processing  experience.  The  chirp  filter 
length  which  produced  the  maximum  signal  peak  was  considered  to  be  the 
"best"  length. 

C.  REGIONALIZATION  OF  EVENTS  BY  CHIRP  FILTER  LENGTH 

Experience  with  signals  recorded  at  ALPA  indicates  that  the 
optimum  length  of  constant  bandwidth  chirp  filters  is  linearly  proportional  to 
epicentral  distance  (Harley,  1971,HeLing  et  al  1972).  It  has  been  found  that  at 
NORSAR  such  a  simple  relationship  between  length  and  distance  does  not 
exist.  Initial  estimates  of  good  chirp  lengths  for  events  from  several  areas 
not  covered  last  year,  e.  g.,  Italy,  Greece,  Iran,  and  Taiwan,  often  were 
quite  wrong,  necessitating  reprocessing  in  these  cases.  Toincrease  processing 
efficiency,  a  table  of  chirp  lenghts  by  region  was  built  up  empirically  and 
refered  to  as  needed.  As  the  number  of  processed  events  increased,  the 


Details  of  the  chirp  filter  derivation  are  outlined  in  Special 

Report  No.  5. 
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table  became  rather  extensive  because  some  regions  required  several  "best" 
lengths.  Thus,  when  a  sufficient  number  of  events  had  been  processed,  maps 
were  produced  showing  the  value  of  the  "best"  chirp  length  in  each  area.  These 
maps  are  shown  in  Figure  V-l  and  V-2  for  the  Love  wave  and  Rayleigh  wave 
chirps  respectively.  The  variat.ons  of  chirp  length  show  distinct  regional 
patterns  which,  as  will  be  discussed  later,  were  the  primary  basis  for  combining 
events  into  groups  of  similar  behavior. 

The  procedure  for  drawing  these  maps  was  as  folio  vs.  The 
optimum  chirp  length  for  each  detected  event  from  the  1971  and  1972  data  was 
plotted  at  the  epicenter  location.*  Contours  were  plotted  to  show  lines  of 
constant  chirp  length.  The  contouring  followed  three  rules: 

•  Chirp  length  should  monotonically  increase  with  distance  along 
any  great  circle  path  away  from  NORSAR. 

•  The  data  were  smoothed  as  to  produce  contours  which  showed 
primary  features  rather  than  each  point-to-point  variation. 

•  Contours  which  were  considered  relatively  reliable  were 
plotted  as  solid  lines.  Contours  based  on  only  a  few  widely 
scattered  events  and  which  were  considered  less  reliable  were 
plotted  as  dashed  lines. 

In  general,  the  maps  for  Love  and  Rayleigh  waves  are  similar, 
as  would  be  expected.  The  data  for  the  Rayleigh  wave  were  much  less  consistent 
than  the  Love  wave  data,  however,  end  hence  were  more  difficult  to  contour. 

For  Love  waves.there  seemed  to  be  very  few  chirp  lengths  which  were  greatly 
di'ierent  from  those  of  nearby  events, while  the  best  chirps  for  Rayleigh  waves 
often  had  lengths  significantly  different.  In  the  one  area  of  the  Caspian  Sea 
and  its  environs,  chirp  lengths  could  not  be  reconciled  with  the  monotonicity 


Chirp  lengths  from  the  1971  events  were  scaled  to  correspond 
to  the  slight  change  in  processing  passbands. 
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OPTIMUM  CHIRP  FILTER  LENGTHS  (IN  SECONDS)  FOR  LOVE  WAVES 
AT  NORSAR  FOR  THE  PASSBAND  0.  025  -  0.  059  Hz 


OPTIMUM  CHIRP  FILTER  LENGTHS  (IN  SECONDS)  FOR  RAYLEIGH  WAVES 
AT  NORSAR  FOR  THE  PASSBAND  0.  025  -  0.  059  Hz 


rule.  Events  from  that  region  consistently  required  longer  chirp  durations 
than  other  events  surrounding  it  but  slightly  further  away. 

It  must  be  stressed  at  this  point  that  these  maps  are  not 
necessarily  unique.  They  are  highly  useful,  however,  for  predicting  what 
chirp  length  is  needed  for  an  event  from  some  particular  location  when  using 
the  0.  025-0.  059  Hz  passband.  They  are  not  unique  for  the  following  reasons: 

•  the  linear  group  velocity  versus  frequency  chirp  only  grossly 
approximates  the  dispersive  characteristics  of  most  events. 

•  the  apparent  "best"  chirp  length  is  strongly  dependent  on  the 
source  energy  spectrum. 

•  the  "best"  length  is  generally  a  nonlinear  function  of  the 
passband  frequencies. 

The  events  used  to  compute  the  SNNR  improvements  of  chirp 
filters  were  grouped  into  13  regions  whose  boundaries  were  suggested  by  the 
behavior  of  the  chirp  length  contours.  These  regions  are  shown  in  Figure 
V-3  and  are  identified  by  a  three  letter  region  descriptor. 

D.  CHIRP  FILTER  SNNR  IMPROVEMENT 

Chirp  filter  SNNR  gains  were  computed  using  187  events  from 
1972  and  45  events  from  1971.  As  mentioned  previously,  the  1971  and  1972 
processing  bands  were  slightly  different,  but  it  was  felt  that  the  SNNR  improve¬ 
ments  would  not  be  significantly  changed. 

The  chirp  filter  SNNR  improvements  are  presented  by  region  in 
Table  V-l.  Included  in  the  table  are  the  lengths  of  the  best  chirp,  the  signal- 
to-noise  ratio  of  the  event  signal  on  the  bandpassed  beam,  the  SNNR  improve¬ 
ments  (in  dB)  and  the  mean  and  standard  deviation  of  the  improvements.  Each 
component  was  treated  independently.  A  non-detection  is  indicated  by  a  dash. 
Gains  were  not  computed  from  all  event-component  detections.  First,  if  none 
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REGION  'ITA'  -  CENTRAL  ITALY,  AUSTRIA,  NCR  THF  R  N  YUGCSLAVIA 


EVENT 

CHIRP  LENGTH 

BANDP ASSET 

SNNP 

SNNR 

IMP 

( 

LOT 

LR  V 

LRR 

LOT 

LR  V 

LFR 

l  OT 

L  °V 

L  po 

AUS*005*04 

HO 

100 

ICO 

1  6.  3 

11.7 

8.9 

0.1 

3  .7 

2.6 

AMS* 169*09 

70 

1  2  C 

120 

12.3 

3  0.4 

14.6 

1.3 

3.2 

1  .9 

I  T  A*035*02 

65 

85 

P  5 

224.4 

62.7 

53.3 

1.9 

1  .  1 

-.3 

IT  A*035*04 

65 

- 

- 

3.8* 

_ 

_ 

0.4 

I T  A*03  5*09 

75 

65 

6* 

131.3 

■*)  .4 

24.5 

0.3 

0.6 

0.’ 

I T  A*03  5*  1  7 

75 

65 

ICS 

114.6 

/9.  5 

25.7 

0."* 

0.8 

-  •  4 

I T  A*035* 1 9 

75 

65 

65 

2R.  9 

H.  5 

4.4  + 

1.4 

0.3 

0  .0 

I T  A*016*05 

65 

95 

55 

22.8 

12.5 

10.7 

0.8 

-.6 

-.3 

IT  A*036*  07 

55 

105 

1  C5 

38.4 

??.■* 

19.8 

0.9 

0.2 

-  .2 

I  T  A*036*  l  5 

75 

95 

P5 

98.3 

34.8 

34.6 

2.0 

0.2 

-1.2 

I T  A*037*0 1 

65 

65 

65 

88.0 

22.1 

24.1 

1.7 

0.3 

-.1 

IT  A*037*2  L 

90 

- 

- 

12.7 

— 

_ 

1  .° 

I T  A*039* 1 2 

75 

105 

105 

,  *3.9 

1  0.4 

10.4 

1  .3 

-.6 

0.3 

ITA*  166* 1 8 

65 

1  10* 

110* 

847.7 

469.4 

381.4 

1.2 

2.8 

*> 

I T A*  166*2  1 

65 

110 

11C 

89.  1 

2  8.7 

16.0 

1  .5 

0.1 

-.3 

YUG* 1 77*  04 

80 

160 

160 

6.3 

15.8 

7.3 

0.9 

3.3 

3  .9 

YUG* 180*01 

80 

175 

150 

42.7 

11.9 

9.P 

1.5 

"»  .0 

1  .8 

MF  AN  SNNR  IMPROVFMFNT  =  1  .  IQ  1.H  0.69 

STANDARD  PFVIATICN  =  0.59  1.52  1.56 

NUMBER  OF  FVFNTS  =  16  14  13 
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REGION  'TUR'  -  GREECE  »Tl  IP  KEY, BLACK  SFA,E  EURf)PF,F  HET>ITFPRANf  ,N 


EVENT 

CHIRP  LENGTH 

PANOPASSET 

SKNR 

SNNP 

IMP 

(OR) 

10T 

LRV 

l  RR 

LOT 

LRV 

LRR 

LOT 

LRV 

LPP 

PL  S/2  63/06 

130 

?  1  7 

217 

0.0 

0.0 

0.0 

1.3 

P.« 

3.3 

Bl  S/263/10 

383 

23  5 

235 

0.0 

0.0 

0.0 

3.2 

3.0 

1  .o 

CAU/2P3/09 

87 

87 

87 

0.0 

0.0 

0.0 

1.0 

0.9 

0.3 

TRS/251/22 

108 

217 

217 

0.0 

0.0 

0.0 

2.6 

2.6 

2.6 

TUR/126/04 

12? 

3  4  F 

34  R 

0.0 

0.0 

0.0 

3.2 

5.3 

4.2 

TUR/143/01 

1CP 

305 

3  C  5 

c.o 

0.0 

0.0 

2.5 

6.5 

6.4 

TUR/161/09 

139 

28^ 

2  83 

c.o 

0.0 

0.0 

4.P 

5.6 

4.9 

PUL*068*22 

*C0 

45C* 

450* 

12.8 

32.1 

14.0 

4.4 

3.6 

2.9 

CAU*079*03 

100 

- 

ICC 

16.6 

- 

5  •  2  * 

-1.5 

- 

1  .3 

C  AU*208* 18 

60 

160 

200 

21.1 

26.  1 

13.2 

0.4 

1.8 

2.0 

CR  E  *01 7*05 

100 

400 

4CC 

6.4 

5.1* 

4.3* 

1  .3 

2.7 

2.9 

CRE*026*12 

- 

250 

- 

- 

2.3* 

- 

- 

-2.1 

- 

CRE*  161 *07 

160 

250 

250 

1  79.0 

161.4 

110.7 

0.5 

3.8 

2.7 

000*020*02 

100 

30C 

3CC 

26. fl 

64.0 

55.5 

1.3 

3.1 

3.9 

GR  E*002*09 

150 

2  PC 

2  PC 

4.7* 

4.1* 

4.9* 

2.0 

2.6 

o.c 

GRE*012*13 

250 

300 

— 

67.5 

126.2 

2.7 

4.1 

_ 

GR  E*03  3*  2 1 

220 

280 

2  PC 

42.  R 

26.0 

10.4 

1.4 

2.0 

2.2 

GRE*0A4*13 

150 

240 

?4C 

43.8 

15.4 

11.3 

1.6 

-.8 

0.7 

GRE*047*00 

150 

280 

240 

51.8 

11.0 

12.5 

2.1 

-.8 

0.2 

GR  E* 157* 1 0 

100 

280 

32C 

138.1 

69.2 

70.5 

-.2 

2.3 

2.3 

GR  E* 167*00 

180 

320 

2P0 

614.1 

562.9 

403.1 

3.1 

5.8 

5.3 

GRE*190*05 

100* 

240 

24C 

11.7 

12.8 

9.5 

0.7 

4.2 

5.5 

GR  E*200* 1 3 

140 

16C 

- 

21.9 

15.1 

- 

3.3 

1.7 

- 

HEC*191*13 

100 

- 

- 

14.1 

- 

- 

2.2 

- 

- 

HE0*199*03 

220* 

?ec 

2  PC 

24.4 

14.8 

18.4 

1.1 

3  .4 

1.2 

HED*205*18 

160 

280 

28C 

9.1 

15.2 

8.7 

1.7 

3.1 

3.8 

SWR*277*08 

50 

40 

6(! 

5.0 

25.6 

13.? 

0.4 

0.3 

0.6 

TRS*160* 12 

100 

280 

320 

13.2, 

5.9* 

5.1* 

0.4 

2.1 

1.6 

TUR*022*17 

125 

250 

25C 

4.4* 

4.6* 

4.0* 

1.9 

3.2 

2.8 

TUR  *  156* 1 6 

125 

250 

3C0 

12.0 

8.4 

7.4 

3.0 

?  •  5 

2.2 

TUR*  170*22 

100 

320 

32  C 

15.2 

8.7 

9.2 

2.8 

2.9 

-1.7 

TUR*173*05 

100* 

200* 

200* 

39.4 

65.6 

3P.5 

0.9 

-.3 

1.1 

TUR*175*04 

120 

240 

2  PC 

31.0 

23.8 

14.0 

3.4 

5.2 

4.0 

TUR  *  186*06 

100 

240 

- 

8.3 

13.7 

- 

2.1 

0.5 

- 

TUR*198*02 

12C 

3  20* 

32C* 

575.4 

335.1 

270.0 

2.3 

3.7 

3.6 
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REGION  'TUR*  -  GREECF  tTURKEYtnlACK  SEA,F  EUROPE, E  MEDITERRANEAN 


EVENT 

CHIRP  LENGTH 

RANOPASSEE 

SNNP 

SNNR 

IMP 

(  DR) 

LOT 

LRV 

LRR 

LOT 

t  RV 

LRP 

LOT 

L  RV 

LRP 

TUR  *206* 10 

100 

?20 

?2C 

34.2 

24.4 

19.2 

1  .9 

3.6 

2.9 

WKZ*233*02 

60 

90 

75 

11.? 

2  7.8 

22.5 

0.1 

0.9 

1.1 

WR  S*204*05 

ac 

120 

12C 

23.4 

8  R.  4 

44.0 

1  .4 

0.5 

1.3 

MEAN  SNNR  IMPROVEMENT  = 

1 .89 

2  .98 

2.52 

STANDARD 

DF VI  ATI  ON  = 

1.37 

1.95 

1.91 

NUMBEP 

EF  EVENTS  = 

3? 

?R 

25 

I] 

; 


l 

j 
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REGION  'CAS 

•  - 

CASPIAN 

SEA  , 

EASTERN 

CALC  A  SOS 

EVENT 

CHIRP  LENGTH 

BANDPASSFD 

SNNR 

SNNR 

IMP 

(  HR) 

LOT 

LRV 

LPR 

LQT 

l  RV 

LRR 

LQT 

1  RV 

LRR 

CAS/1  35/04 

287 

65 

65 

0.0 

0.0 

0.0 

3.6 

1.2 

1.0 

C  AU/ 262/06 

- 

392 

- 

- 

0.  0* 

— 

— 

0.0 

— 

CAU/288/17 

305 

A  5  3 

653 

0.0 

0.0 

0.0 

2.5 

2.3 

3  .4 

WKZ/356/06 

261 

1  76 

1  76 

0.0 

0.0 

0.0 

1.9 

1.6 

1.2 

CAU*166*00 

350 

300 

300 

22.6 

76.6 

66.0 

0.2 

1.3 

1.0 

I R  A*0 1 8*2 1 

30C* 

60C 

6C0 

63.1 

16.  7 

13.3 

-1.2 

3.9 

1.5 

MEAN  SNNR  IMPROVEMENT  = 

2.00 

2.03 

I  .63 

STANDARD 

DEVIATION  = 

1.34 

1.13 

1.01 

NUMBER 

OF  EVENTS  = 

6 

5 

5 
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REGION  'IRA*  -  SOUTHERN  IRAN 


EVENT 

CHIRP  LENGTH 

BANDPA  SSF  0 

SNNR 

SNNR 

f  VP 

I  OP) 

LOT 

LRV 

LPR 

LOT 

LRV 

LPR 

LOT 

l  PV 

LRR 

IR  A*006*  09 

3  00* 

420 

_ 

47.4 

31.1 

_ 

-2.4 

0.8 

_ 

IRA*014*22 

25C* 

56C 

560 

93.6 

33.4 

19.6 

-2.1 

-.8 

?.l 

IRA*041*09 

200 

440 

440 

21.8 

7.7 

6.6 

-.8 

1.7 

2.6 

IR A*041*16 

450 

400 

4CC 

11.8 

23.5 

8.5 

0.9 

0.0 

1.0 

I R  A*068*2 1 

350 

- 

400 

19.4 

- 

cc 

• 

rv 

cv 

0.2 

- 

2.7 

IR  A*1 55  +  08 

160* 

360 

— 

11.3 

3.9* 

0.4 

1.7 

IRA+156+03 

- 

540 

540 

- 

3.3* 

4.7* 

- 

5.5 

0.6 

IRA*157*1 1 

160* 

300 

360 

45.4 

8.  8 

8.R 

-2.3 

1.7 

-1.1 

IR  A* 162* 19 

160* 

48C 

480 

63.6 

20.7 

15.9 

-1  .6 

2.3 

3.5 

IRA* 164*13 

100 

350 

35C 

367.0 

68.4 

76.3 

0.6 

2.3 

1.5 

IRA*166*00 

160 

540 

480 

444.9  ' 

14  8.0 

120.0 

-.4 

3.4 

3.5 

IRA*  168*23 

160* 

300 

300 

22.2 

24.4 

8.9 

-2.4 

0.7 

1  .2 

I RA*175*0fl 

160* 

30C 

3C0 

100.4 

25.3 

21.3 

-.8 

3.8 

2.9 

IR  A*  104* 1 2 

200* 

- 

- 

1772.2 

- 

- 

-1.6 

— 

— 

IR  A* 1 84* 14 

90 

350 

3CC 

2.9* 

2.4* 

2.6* 

1.7 

2.2 

0.2 

IR A*185*02 

— 

300 

30C 

— 

32.5 

21  .5 

3.2 

3.4 

I R  i* 185*12 

- 

350 

- 

- 

7.7 

- 

— 

0.5 

— 

IP  A*  18  5*2  1 

100 

350 

3CC 

242.6 

97.  8 

93.3 

i .  3 

4.2 

3.9 

!RA*187*16 

100 

250 

26C 

15.6 

2  1  .R 

21  .4 

3.8 

2.4 

1  .3 

IRA*187*21 

130 

— 

- 

3.2* 

- 

- 

1.4 

- 

- 

IRA*188*05 

— 

350 

350 

6.1 

5.7* 

1.6 

0.7 

IR A*193*22 

- 

200 

200 

- 

54.2 

57.0 

— 

3.3 

2.8 

IRA*196* 13 

100 

400 

- 

21.3 

7.3 

— 

0.9 

0.6 

— 

I RA*196*17 

100 

- 

- 

6.1 

- 

- 

1.1 

— 

— 

IRQ*060*08 

250 

4  0C* 

4  CO* 

11.7 

3.5 

4.4* 

1.1 

1.7 

-3  .4 

MFAN  SNNR  IMPROVEMENT  = 

0.87 

1  .86 

2.23 

STANDARD 

DEVIATION  = 

1.25 

1.42 

1.35 

NUMBER 

OF  F  VENTS  = 

10 

17 

14 

V- 13 
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REGION  'PAK*  -  PAKISTAN,  EASTERN  IRAN 


EVENT 

CHIRP  LENGTH 

BANDPASSEP 

SNNR 

SNNR 

IMP 

(08) 

LQT 

I.R  V 

LRR 

LOT 

LRV 

LRR 

LOT 

L  RV 

LRR 

IRA*029*09 

200 

360 

36C 

4.2* 

6.1 

3.9* 

-1.0 

-1.7 

0.9 

IR  A* 193*  15 

220 

- 

- 

4.2* 

- 

- 

2.9 

- 

- 

I WP*02B* 10 

240 

250* 

2  5  C* 

3  8.9 

31.3 

24.4 

0.0 

-.3 

-2.1 

P  AK  ♦  1 57*  1 1 

240 

1 60 

320 

26.  1 

14.8 

14.0 

3.2 

1.1 

1.1 

P  AK  *  1 62*  1 1 

— 

300 

25C 

- 

46.6 

32.6 

- 

2.3 

1.4 

P  AK* 1 79* 1 0 

280 

240 

16C 

48.5 

132.  1 

IT 

• 

m 

CD 

0.7 

2.1 

1.7 

PAK*  195*18 

- 

240 

- 

- 

6.5 

- 

- 

2.8 

- 

MEAN  SNNR  IMPROVEMENT  =  1.30  1.41  U41 

STANDARD  DEVIATION  =  1.70  1.61  0.28 

NUMBER  OF  EVFNTS  =  353 
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REGION  'TO? 

•  - 

TADZHIK 

,  KIRGI7 

,  HINDU  KUSH, 

N  A  EG  AM  STAN 

EVENT 

CHIRP  LFNGTH 

FANDPASSEO 

SNNR 

SNNR 

IMP 

<PP) 

LOT 

LRV 

LPR 

LOT 

LRV 

LRR 

LCT 

LRV 

LPo 

KRG/301 / 13 

122 

383 

383 

0.0 

0.0 

o 

. 

o 

l  .5 

2.9 

1.8 

TCZ/147/00 

348 

292 

- 

0.  0 

0.0 

- 

2.1 

3.6 

— 

TDZ/274/16 

348 

479 

479 

0.0 

0.0 

0.0 

-3.2 

3.7 

3.8 

AFG*059* 18 

250 

- 

400 

6.6 

- 

13.5 

3.8 

- 

4  .  ^ 

AFG* 181*03 

300* 

40  C* 

4C0* 

54.1 

66.6 

50.2 

-.6 

-1.0 

0.9 

HNK*053*U8 

250 

35C 

25C 

7.8 

4.5* 

3.7* 

2.1 

0.° 

2.0 

HNK* 177*07 

380 

- 

- 

48.7 

- 

- 

-1  .4 

- 

— 

KRG*006*06 

80 

100 

130 

4.1* 

9.7 

6.7 

-1.6 

0.3 

7.2 

T  AD*077*09 

400 

350 

3  50 

69.3 

136.3 

121.0 

1.3 

1.0 

1.7 

TDZ*005*12 

450 

450 

450 

6.7 

6.6 

6.6 

2.3 

O.Q 

0.5 

HNK*  178*20 

_ 

700 

70C 

- 

4.6* 

5.0* 

- 

3.2 

2.6 

HNK* 179* 1 5 

- 

200 

2CC 

— 

107.7 

76.2 

2.1 

2.5 

MEAN 

SNNR  IMPROVEMENT  = 

1  ,07 

2.07 

2.38 

STANDARD 

DF VI  ATI  ON  = 

2.25 

1.38 

1.26 

NUMRFR 

DF  EVENTS  = 

8 

7 

7 
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REGION  'TIB'  -  S  SINKIANG,  TIBET,  N  BURMA,  HIMLAYAN  MOUNTAINS 


EVENT 

CHTRP  LFNGTH 

PANDPASSFD 

SNNP 

SNNP 

IMP 

(DP) 

LOT 

LRV 

LPP 

LOT 

LRV 

LRR 

LOT 

LRV 

LRR 

SIN/219/15 

17ft 

305 

3C5 

o 

. 

o 

0.0 

O 

. 

o 

3.1 

3.0 

0.3 

S IN/241/ 1 5 

- 

17ft 

- 

- 

0.0* 

- 

- 

0.7 

- 

TIB/123/00 

196 

305 

3C5 

C.O 

0.0 

0.0 

4.5 

1.2 

3.5 

TIB/155/20 

305 

ft  79 

ft7S 

0.0 

0.0 

0.0 

3.2 

1.0 

1.6 

TIB/302/17 

217 

87 

e7 

0.0 

0.0 

0.0 

2.3 

1.0 

1.2 

RIIR  *  160*  16 

350 

ftOO 

7C0 

19.7 

26.5 

13.7 

3.2 

-1.4 

-.4 

CH I  *  1 5ft*  1 6 

250 

ft20 

ft20 

42.8 

16.5 

13.8 

2.4 

3.8 

3.1 

IHO* 154*20 

250 

ftftC 

- 

10.8 

7.5 

- 

1.9 

-.7 

- 

T I B*075*06 

300 

300 

30C 

ft  0.  7 

37.1 

24.0 

0.7 

2.0 

2.R 

T 1 R* 160*  23 

120 

100 

175 

16.4 

5.R* 

5.2* 

l.o 

1 .0 

2.3 

TIB*  170*04 

2ft0 

3  60 

36C 

35.2 

36.2 

14.3 

-.8 

3.2 

1.1 

T I B* 195*05 

200 

- 

- 

6.7 

- 

- 

3.2 

- 

— 

TIP*19fl*02 

250 

- 

• 

66.4 

- 

- 

3.4 

— 

— 

T I P*193*03 

200 

- 

- 

10.6 

* 

- 

2.0 

— 

- 

T I B*20ft* 16 

150 

200 

2  CO 

5416.3 

936.4 

511. R 

2.R 

1.3 

-1.0 

T I B*  204*2 1 

150 

— 

— 

27.9 

- 

- 

3.7 

- 

- 

T  I  P*205*23 

150 

- 

- 

R.  ft 

- 

- 

3.5 

- 

— 

T I R*  206* 1 ft 

150* 

200* 

200* 

20.9 

20.7 

15.3 

-.7 

1.2 

O.R 

MF AN  SNNR  IMPROVEMENT  =  2.56  i.43  1.35 

STANDARD  DEVIATION  =  1.27  l.ft5  1.5B 

NUMBER  CF  EVENTS  =  1ft  ID  9 
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RFGION  *CH  I  *  -  PENTRAL  CHINA 


EVENT 

CH  1 R 

°  LENGTH 

B6  NDP ASSC  D 

SNNR 

9NNR 

I  MP 

(  OR) 

LOT 

LRV 

LPR 

LOT 

LRV 

L  R  P 

LOT 

1  RV 

L  RR 

CHI/229/09 

152 

6  79 

479 

0.0 

0.0 

0.0 

2.3 

4.0 

4.3 

CHI/229/17 

152 

196 

1  c6 

0.  0 

0.0 

0.0 

1.6 

1.0 

1  .4 

CHI/258/07 

1 76 

471 

471 

0.0* 

0.0* 

0.0* 

2.5 

3.2 

<.l 

01*036*07 

180 

550 

550 

37.7 

1  ->.3 

Q.  7 

2.3 

1  .Q 

?.3 

CHI*  189*2  3 

380 

- 

- 

15. 9 

- 

1  .5 

CH I *203*  16 

310 

40C 

4  CO 

3  P.  3 

69.6 

29.9 

5.8 

3.6 

4.6 

MON*  153*1 1 

200 

- 

- 

9.  9 

— 

_ 

1  .7 

YUN*05  7* 1 8 

300 

30C 

fcCO 

45.7 

3^.2 

20.3 

0.5 

-.3 

-1.4 

MEAN  SNNP  I M PRTVF  ME  NT  =  2 . 24  2.22  2.26 

STANDARD  DEVIATION  i  1.67  1.71  ?.43 
NUMBER  flF  EVFNTS  =  7^6 
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REGION  'SIN*  -  SOUTHWESTERN,  WESTERN,  NOPTHWFSTf RN  SINKIANG 


EVENT 

CHIRP  LENGTH 

BANDPASSEC 

SNNP 

SNNR 

IMP 

1  OP ) 

LOT 

LRV 

LPR 

LOT 

LRV 

LPR 

LOT 

LRV 

LRR 

SIN/166/22 

87 

174 

174 

0.0 

0.0 

0.0 

1.6 

1.2 

0.0 

S IN/166/23 

95 

217 

217 

0.0 

0.0 

0.0 

1.4 

-.6 

-•? 

S IN/170/17 

ICE 

217 

217 

0.0 

0.0 

0.0 

0.7 

1.3 

1  .9 

S IN/221/01 

104 

78 

78 

0.0 

0.0 

0.0 

1.9 

0.6 

0.3 

CRS/236/16 

113 

235 

- 

0.0 

0.0 

- 

2.6 

2.5 

- 

SIN/273/12 

104 

174 

1  74 

0.0 

0.0 

0.0 

1.7 

1.5 

0.4 

RA I *058*22 

100 

300 

3C0 

7.7 

12.2 

8.4 

1  .2 

2.7 

2.5 

EKZ  *078*07 

- 

5C 

50 

- 

4.7* 

3.3* 

- 

0.4 

1.3 

KR  G*028*20 

80 

32C* 

32  C* 

9.7 

4.5 

3.7* 

0.8 

-1.8 

-1.3 

SIN*002*10 

80 

100 

ICC 

70.4 

75.8 

62.4 

2.5 

1.1 

1.0 

S I N*042*05 

ICO* 

100 

100 

53.3 

70.6 

52.6 

0.9 

l.l 

0.4 

SIN*047*23 

80 

150 

150 

17.1 

31.8 

22.4 

0.8 

0.5 

0.8 

S IN*064*04 

100 

200 

20C 

15.3 

44.6 

19.1 

-.2 

0.7 

1.0 

S  IN*  154*06 

80 

175 

- 

29.0 

5.8* 

- 

1.5 

1.2 

- 

S IN*  187*01 

60 

175* 

175* 

256.0 

78.2 

53.0 

1.1 

1.5 

2.3 

S IN*  1 87*04 

140 

175 

175 

8.0 

12.7 

8.4 

4.2 

3.2 

2.5 

SIN*  192*19 

160 

175 

175 

6.9 

35.6 

26.4 

3.2 

1.6 

1.1 

SIN*200*03 

100 

- 

- 

7.6 

- 

- 

0.7 

- 

- 

MEAN  SNNR  IMPROVEMENT  =  1.60  1.3*-  0.97 

STANDARD  DEVIATION  =  1.08  1.01  0.80 

NUMBER  OF  EVENTS  *  IS  13  12 
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REGION  '  EK  Z  '  -  EASTERN  KAZAKH  TEST  AREA,  CENTRAL  KAZAKH 


EVENT 

CHIRP  LENGTH 

BANOPASSFC 

SNNR 

SNNR 

IMP 

(  DP) 

LOT 

LRV 

LPR 

LOT 

LRV 

LPR 

LCT 

LRV 

LPR 

EKZ/145/04 

- 

113 

0.0* 

4-0 

FKZ/157/04 

91 

200 

20C 

0.0 

0.0 

0.0 

1 . 3 

1 . 8 

0.2 

EK  * /1 70/04 

52 

139 

135 

0.0 

0.0 

0.0 

0.5 

0.3 

1  -  P 

EKZ /282/06 

- 

1  35 

- 

— 

0.0* 

2-7 

EKZ/294/06 

39 

122 

122 

0.0 

0.  0 

0.0 

0.9 

O.R 

-1.6 

LKZ/333/C6 

52 

65 

69 

0.0 

0.0 

0.0 

1.4 

0.5 

1  .  1 

EKZ/364/06 

65 

P  7 

P7 

0.0 

0.0 

0.0 

1  .3 

1  .  1 

EKZ  *070*04 

50 

100 

ICO 

10.0 

5.9 

6.6 

1  .0 

1-4 

1  .  r 

1  ? 

EKZ  *088*04 

130 

50 

50 

2.  8* 

4.0* 

3.2* 

4.2 

0.1 

1  •  < 

f)  -  1 

EKZ  *229*0  3 

• 

60 

eo 

- 

6.3 

5.4* 

0.? 

0.1 

EKZ*  307*0 1 

- 

75 

29.  0 

1  ft 

CKZ*345*04 

50 

50 

75 

77.3 

34.0 

29.  1 

o .  q 

1  1 

K A7  *  18 1 *00 

— 

— 

ICO 

- 

^.7* 

•  * 

V*  •  ’ 

1  •  l 

1.4 

MEAN  SNNR  I MPROVFMFNT  = 

0.87  0 

.97 

0  .71 

STANDARD 

DEVIATION  * 

0.61  0 

.60 

1.12 

NUMBER 

OF  EVENTS  * 

7 

9 

7 

V- 19 
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REGION  ' T A I •  -  TAIWAN,  S  RYUKYU  ISLANOS,  EAST  CHINA  SFA 


EVENT 

CHIRP  LENGTH 

bandpassfo 

SNNR 

SNNR 

IMP 

(  OR) 

LOT 

LRV 

LRR 

LOT 

LRV 

LRR 

LOT 

LRV 

LRR 

CHI*028*04 

60C* 

700 

700 

11.3 

7.5 

6.9 

-2  .4 

0.0 

•>  ? 

RYU*  155*02 

300 

50C 

500 

108.3 

101.7 

66.5 

1 .2 

3.6 

w  •  < 

RYU* 197*02 

24C 

52C 

520 

24.2 

48.0 

34.1 

3.1 

3.0 

9  .  c  » 

RYU*209* 16 

36C 

64C 

660 

261.4 

106.8 

68.5 

2.9 

3.9 

’  •  o 

T  A  I  *006*05 

500* 

1000 

- 

5.0 

3.2* 

2.9 

1.8 

1  •  ; 

T  A  1*006*  1 2 

600* 

60  C 

60C 

97.9 

64.0 

50.2 

2.4 

1 . 5 

'l  1 

T  A  1*006*06 

500* 

600 

600 

19.7 

10.5 

13.5 

0.3 

0 . 1 

.  1  •  l 

•  o 

TAI*010*05 

700 

800 

PCO 

50.6 

20.0 

17.7 

0.2 

4.0 

a  q 

T  A  I* 1 78*08 

300 

5CC 

5CC 

1  58.7 

c4 . 8 

81.4 

2.3 

4  2 

t  «  n 

U 

TAI*182*18 

•>oo* 

5C0 

500 

548.  8 

153.3 

80.  1 

3.7 

“it 

0.8 

•  j 

2.4 

T  A l ♦ 19  5*  2  3 

3C0 

_ 

20.1 

a  n 

TAI*198*13 

300 

600 

600 

30.0 

55.8 

26.4 

3.6 

3.8 

4.1 

MEAN  SNNR  IMPROVEMENT  = 

2.74  2 

.59 

2.89 

STANDAPO 

OF  VI  ATI  ON  = 

1.85  1 

.81 

1.75 

NUMRER 

OF  EVFNTS  = 

7 

10 

10 

V-20 
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REGION  'KUR'  -  KURILE  ISLANDS 


EVENT  CHIRP  LENGTH 


LOT 

LRV 

LRR 

KUR/135/21 

828 

°06 

8C6 

KUR/146/01 

566 

762 

762 

KUR/147/16 

479 

828 

828 

KUR / 152/21 

566 

697 

697 

KUR / 190/ 1 6 

828 

1C02 

1CC2 

KUR/191/03 

741 

958 

95e 

KUR/191/09 

523 

871 

871 

KUR  /  2 13/02 

459 

893 

893 

K  AM*078* 18 

980 

- 

— 

KUR  *00 1*16 

700 

- 

• 

KUR*00 1  *  1 8 

700 

KUR+005*02 

700 

O0C 

9C0 

KUR*009* 14 

- 

- 

9C0 

KUR*022*01 

- 

700 

1150 

KUR*028*23 

980 

- 

- 

KUR  *046* 1 6 

860 

8  50 

850 

KUR*049* 1 8 

860 

850 

850 

KUR*054*03 

620 

1000 

1000 

KUP  *054*2  7 

740 

100C 

10CC 

<kUR*054*37 

62  C 

1000 

10C0 

KUR*055* 10 

980 

100C 

ICCO 

KUR*056*22 

620 

700 

700 

KUR*057*02 

980 

700 

700 

KIJR*057*05 

980 

850 

850 

KUR*063*23 

740 

1000 

10C0 

KUR*070*  06 

620 

850 

KUR*077*  07 

980 

- 

- 

KUR*153*00 

500* 

700 

7CC 

KUR  *  169*  1 9 

740 

1000 

1  300 

KUP  *  171*18 

9eo 

700 

11  50 

KUR  *171*22 

— 

1C00 

KUR ♦  190*2  1 

500 

650 

10C0 

KUR*193*06 

- 

850 

- 

K'JP  *  194*00 

860 

05C 

ICCO 

KUR* 195* 1 5 

500 

700 

7C0 

FROM  JAP**N  TC  KAMCHATKA 


RANDPASSED 

SNNP 

SNNP 

IMP 

(  OP) 

LOT 

l  R  V 

LRR 

LOT 

LRV 

LRR 

0.0 

0.0 

0.0 

-.6 

4.8 

4.1 

0.  0 

0.0 

0.0 

5.7 

6.8 

6.2 

0.0 

0.0 

0.0 

6.3 

4.8 

7.1 

0.0 

0.0 

0.0 

2.7 

5.9 

4  .9 

0.  0 

0.0 

0.0 

6.9 

4.4 

4  .6 

0.0 

0.0 

0.0 

6.2 

4.0 

7.2 

0.0* 

0.0* 

0.0* 

0.5 

O.P 

0.5 

0.0 

0.0 

0.0 

4.7 

P.7 

7.9 

3.6* 

- 

- 

1  .8 

- 

— 

3.4* 

— 

— 

1.2 

- 

- 

- 

4.5* 

— 

0.0 

4.7* 

4.7* 

3.6* 

2.6 

6.1 

6.2 

- 

- 

2.2* 

- 

- 

2.0 

- 

3.5* 

4.4* 

- 

1.0 

0.0 

3.4* 

— 

- 

-4.3 

- 

- 

3.4* 

4.  8* 

2.9* 

0.5 

5.0 

4.2 

4.2* 

4.9* 

3.8* 

3.8 

-.8 

1.7 

36.8 

26.6 

23.3 

2.9 

0.7 

0.7 

25.0 

19.6 

16.9 

3.8 

2.2 

2.7 

30.7 

24.4 

21  .  7 

3.1 

1.1 

0.9 

^6.7 

134.4 

90.4 

3.5 

5.7 

5.7 

3.5* 

4.0* 

3.7* 

1.8  ■ 

-1.5 

0.7 

46.6 

132.5 

102.1 

2.7 

1.8 

0.6 

6.8 

7.4 

6.6 

3.4 

7.1 

7.0 

4.5* 

1  1.  P 

6.4 

3.0 

6.1 

10 .0 

’.5* 

3.0* 

— 

2.3 

l.P 

11.4 

- 

- 

4.6 

- 

- 

8.5 

6.0* 

6.6 

3.6 

4.2 

5.2 

11.2 

10.7 

5.7* 

2.5 

3.7 

0 

14.3 

20.7 

9.  ? 

0.2 

0.6 

7 

- 

5.8* 

_ 

• 

3.7 

— , 

2.7* 

5.0* 

3.0* 

3.9 

1.6 

4.5 

- 

51.1 

- 

- 

7.0 

- 

23.3 

53.5 

3°.  4 

4.5 

7.7 

7.4 

4. 8* 

9.9 

6.5 

4.7 

4.7 

2.7 
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REGION  *KUR '  -  KURILF  ISLANOS  p ROM  JAPAN  TC  KAMCHATKA 


EVENT 

CHIRP  IFNGTH 

BANDPA  SSEO 

SNNR 

SNNR 

IMP 

(DB) 

LQT 

LRV 

LPR 

LOT 

LRV 

LRR 

LQT 

LRV 

LPR 

KUR  *  2 1 1*  2 1 

50C 

700 

7CC 

55.6 

81.1 

34.6 

2.4 

1.7 

2.8 

MEAN  SNNR  IMPROVE 

ME  NT  * 

3.64 

4.48 

4.77 

STANDARD 

DEVIATION  * 

1  .98 

2.46 

2.84 

NUMRFR 

OF  EVENTS  * 

18 

?0 

19 
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REGION  ‘KAN*  -  KAMCHATKA  PF N 1  NSL'LA  ,  KCMANOCR  SK Y  ISLANDS 


EVENT 

CHIRP  LENGTH 

RANOPASSFO 

SNNR 

SNNR 

1  Vp 

inp) 

LOT 

LRV 

L  PR 

LOT 

LRV 

LRP 

LOT 

LRV 

LRR 

KAM/166/14 

610 

1046 

1C46 

0.0 

0.  0 

0.0 

-.1 

1.3 

O.R 

KOM/ 148/10 

522 

697 

697 

0.0 

0.0 

0.0 

0.0 

4.6 

4.3 

K  AM*003*06 

660 

700 

70C 

12.5 

1  6.6 

14.3 

5.3 

2.0 

3.6 

K  AM*003* 19 

700 

- 

- 

3.5* 

- 

- 

0.8 

- 

- 

K  AM*004*02 

600 

OOO 

eco 

3.6* 

3.0* 

2.7* 

-1.2 

2.3 

2.1 

KAM*005*16 

620 

e5c 

P5C 

3.6* 

■'.4* 

3.7* 

3.7 

1.8 

-.3 

KAM*009*14 

620 

- 

- 

3.7* 

- 

- 

-1.1 

- 

- 

K  AM *0 1 2* 20 

5CC 

1050 

l’CO 

16.1 

13.3 

11  .? 

2.4 

1.4 

*  •  6 

KAM*042*2l 

740 

1000 

850 

6.6 

6.2 

6.9 

-.2 

2.0 

-1.6 

K  AM*052*22 

- 

1000 

— 

— 

3.9* 

— 

- 

2  .q 

— 

K  AM  *  1 56*  07 

— 

1150 

- 

- 

5.3* 

- 

- 

1.4 

- 

K AM*157*04 

620 

1300 

ICC 

8.2 

R  •  3 

6.1 

2.6 

0.5 

1  •  M 

K  AM* 163* 14 

- 

10  30 

- 

- 

5.6* 

- 

- 

0.3 

- 

K  AM* 1 77* 17 

620 

IHO 

- 

6.1 

5.6* 

- 

3.* 

4.1 

- 

K  AM* 180*14 

— 

I’OC 

100C 

— 

3.3* 

6.0 

— 

1.8 

0.4 

K  AM*  186* 1 3 

— 

P  50 

— 

- 

4.0* 

- 

- 

4.4 

- 

K  AM* 189*05 

- 

noo 

- 

- 

3.2* 

- 

- 

5.q 

- 

KAM*193*0R 

50C 

- 

- 

4.0* 

- 

- 

2.7 

- 

- 

K  AM* 19 7* 1 3 

- 

1000 

- 

- 

3.2* 

- 

- 

1.3 

- 

K  AM*  206* 1 3 

- 

70C 

7CC 

— 

8.2 

6.4 

— 

2.2 

2.3 

K0M*044*22 

74  0 

700 

700 

3.3* 

3.7* 

3.2> 

3.0 

2.? 

2.3 

KOM*  15  3*21 

960 

70C 

700 

6.3 

5.9* 

6.1 

-1.4 

2.7 

->  .f 

KUR*209*00 

500 

1150 

1150 

104.4 

68.1 

68.  R 

6.1 

1  .9 

’  .4 

VF  A  N  SNNR  I M  TROVE  ME  NT  =  1.9q  1  .9°  ?.A? 

STANDARD  prviATION  =  2.6?  l.?0  2.25 

NIJHOtR  rp  EVENTS  =  9  8  10 
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CHIPP  MAT  rHcn  Pll.TTp  I  MPPPVr  wr  NT  IN  *  1  G  N  A  l  t\n  ISr  RATIO 

°  a r«  f  i  f,  nr  i  h 

RFCItlN  'nTHTP  •  -  0  T  H  r  w  r\/rMS  NCT  r,wr  |_prg 

FVFNT  CHIRP  IfUr.TH  PAN"PASSFr  SN'NP  S\NR  jwp  (i>h) 

LOT  I.RV  Lec  LOT  L°V  I  RP  LOT  \  «V  I  PR 

CFI/16*/10  163  ’Al'  ^Ay  0.0  0.0  0.0  O.a  ?.?  ->.y 

CMI/P49/21  30?  7«A  7fc4  0.0  0.0*  0.0*  2.1  P.7  } . o 

FRS/165/13  17A  ’R.?  OR?  C,.P  0.0  O.Q  4.7  1.0  1.1 

FRS/266/21  152  305  3C5  0.0 *  0.0  0.0*  2.8  -.6  0.4 

RKl  *0',‘>»03  ?0C  30C  OCC  0.7  3.0  3.6*  2.1  1.1  0.? 

CFI*030*03  150  ^OC  ?6 0  16. «  ?7.o  9.R  p.O  J.6  o.O 

CPS*015*18  ?0C  OCO  ‘‘CC  2P.6  10.6  l?.l  '^.o  ^  g*p 

FRS*172*CO  ACC  -  -  y.7  -  _  _.7 

K  AM*05 1*20  -  0  6  C  F^C  .  2.0  3.1*  -  ?.o 

1  CJM  ♦  1 5  7*  1  ??0  1 6  C  ICC  lb. 2  33.3  20.,?  P.P  ?.4  ?!b 

NRS*058* 10  2?0  ’5C*  2CC*  41.0  0?.0  6O.7  *.o  0.4  j  .4 

NRS*058*17  22C  '’CO  ’CO  A.^*  7.4  (■,.?  p.I  0.6  1.6 

RYU*  106*18  -  *80  pF0  .  3.2  3.-*  _  ?#4  ,J#4 

RIP *013*1 7  ACC  QC0  FCC  167. P  0,,,.i  PP^.3  1 .A  6.6  A. A 

S I P*0 1A*03  800  600  6C0  3.?  5.1  6.0*  3.6  A.l  4.8 

U  AR  *  1 PO*  OO  300  A  6C  A60  1  66. 4  73.0  90.3  -.0  6.Q  A. 7 

WF  A  N  lyp^nvrPFM  =  2.0A  2.?6  ?m?/. 

STANDARD  f'c  V I  AT!  CN  -  1  .65  1  .06  1.56 

NllMH L  R  PF  E  VF  NTS  =  1?  10  0 
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of  the  chirps  applied  to  a  particular  component  was  considered  to  be  the 
optimum  length,  an  asterisk  was  placed  next  to  that  length.  Second,  if  the 
(amplitude)  SNNR  of  the  bandpassed  signal  was  less  than  6.  0  (15.  6dB)  it  was 
considered  to  be  too  noisy  to  give  unbiased  results  and  an  asterisk  was  placed 
next  to  the  SNNt  value.  *  If  an  asterisk  of  either  type  was  placed  on  that 
event-component  entry,  its  improvement  was  not  used  to  compute  the  average 
improvement. 

The  SNNR  improvements  obtained  by  chirp  filters  range  from 
about  0.8  dD  for  region  EKZ  to  about  3.  5-4.  5  dB  for  region  KUR.  Median 
improvement  is  2  dB.  Overall,  chirps  applied  to  Love  waves  yield  slightly 
less  gain  than  chirps  applied  to  Rayleigh  waves  (vertical).  There  are  signifi¬ 
cant  differences  however  in  LQ- LR  behavior  from  region  to  region.  For 
example,  in  KUR,  TDZ,  and  IRA,  LQ  gains  were  1  dB  smaller  than  LR  gains; 
while  in  TIB  the  reverse  was  true.  Five  regions  had  LQ>  LRV  ga  ns  and 
seven  regions  with  LQ<  LRV  gains. 

Physical  considerations  indicate  that  the  chirp  gains  for  the 
vertical  and  radial  components  should  be  similar;  hence  the  ratio  of  vertical 
to  radial  gains  should  be  near  unity.  The  average  ratio,  computed  over  145 
events,  is  1.009  (0.04  dB). 

The  standard  deviations  (sigmas)  of  the  SNNR  improvements 
are  in  many  cases  larger  than  the  mean  improvements,  although  most  values 
are  about  0.  6  of  the  mean.  Large  means  tend  to  have  large  sigmas  and  vice- 
versa.  Some  regions  and  wave  components  which  show  relatively  small  sigmas 
include  ITA  (LQ),  TDZ  (LR),  and  EKZ  (LQ,  LRV).  Regions  with  a  relatively 
large  sigma  include  ITA  (LR),  IRA  (LQ),  TDZ  (LQ),  TIB  (LR),  and  KAM  (LQ). 


Bandpassed  SNNR's  of  th^  1971  events  are  shown  as  0.  0 
however  they  also  were  checked  for  actual  S  NNR. 
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The  ratio  of  the  standard  deviations  of  the  radial  and  vertical 
components  of  LR  should  also  be  approximately  equal  to  unity.  This  is 
generally  the  case  with  nine  regions  having  values  between  0.  85  and  1. 15. 

Three  regions,  PAK,  EKZ,  and  KAM  showed  unusually  large  or  small  sigma 
ratios . 

The  sole  purpose  of  using  matched  filters  is  to  lower  the 
detection  threshold  so  that  more  smaller  events  can  be  detected.  The  surface 
wave  magnitudes  assigned  to  such  events  include  a  correction  factor  to  account 
for  the  matched  filter  gain.  What  factor  is  used  should  be  a  funcfion  of  mean 
gain  and  the  standard  deviation  computed  so  as  to  satisfy  whatever  probabilistic 
detection  model  is  used.  In  this  report,  magnitude  correction  factors  were 
computed  using  the  mean  gain  of  each  region. 

The  regionalization  scheme  based  on  chirp  length  does  not 
reflect  as  clearly  any  SNNR  gain  differences  between  regions,  primarily 
because  of  the  large  gain  variances.  For  example,  the  small  region  TDZ,  which 
had  long  chirp  lengths,  has  small  LQ  gains  and  larger  LR  gains.  On  the 
other  hand,  the  three  regions  contiguous  to  TDZ,  PAK.  TIB,  and  SIN,  all  have 
large  LQ  gains  and  smaller  or  equal  LR  gains.  The  event  KRG*006*06  in 
region  TDZ  actually  lies  on  the  TDZ-SIN  regional  boundary.  It  has  the  short 
chirp  lengths  appropriate  to  region  SIN  even  though  its  SNNR  gain  behavior,  con¬ 
sidering  its  low  SNNR,  is  more  like  TDZ  than  SIN.  A  second  example  is  the 
region  KAM,  which  shows  a  distinctly  different  mean  gain  and  standard 
deviation  than  the  main  Kurile  Island  region  KUR  even  though  these  two  regions 
have  similar  chirp  lengths. 

The  region  EKZ  was  differentiated  from  SIN  primarily  because 
events  there  were  mostly  presumed  explosions  rather  than  because  of  inherent 
differences  in  chirp  behavior.  EKZ  events  have  the  smallest  variance  in  gain 
for  LRV  (LQ  gain  variance  was  the  second  smallest)  of  any  region  observed. 

The  presumed  source  mechanism  and  small  geographical  distribution  are 
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most  likely  the  reasons  for  this.  Similarity  of  source  mechanism  and  location 
also  contribute  to  the  small  LQ  gain  variance  for  the  ITA  events.  If  only  the 
Italian  events  are  considered  in  this  region,  the  LR  gains  would  also  have  less 
variance.  Interestingly,  the  Austrian  and  northwestern  Yugoslovian  events 
show  similar  LQ  gains  but  distinctly  larger  LR  gains  than  do  the  Italian  events. 

Chirp  gains  obtained  at  NORSAR  are  similar  to  those  obtained 
at  ALP  A  (Strauss,  1973).  Typical  gains  forevents  at  comparable 
distances  were  generally  fairly  close,  except  for  a  few  areas  where  most  gain 
differences  could  be  attributed  to  differences  in  chirp  length.  For  example, 

ALPA  showed  larger  gains  both  in  eastern  Kazakh  and  in  Taiwan.  Best  (Rayleigh) 
chirp  lengths  for  Kazakh  events  was  about  700  seconds  for  ALPA  but  only  50-100 
seconds  for  NORSAR.  For  Taiwan,  ALPA  lengths  were  about  800  seconds 
compared  to  300  seconds  for  NORSAR. 

The  application  of  chirp  filters  to  the  data  produced  eleven 
event  detections  wh^re  the  surface  waves  had  been  undetected  on  the  bandpassed 
beam.  Ten  of  these  events  had  m^'s  between  3.  6  and  4.  2,  which  brackets 
the  50%  detection  level*,  and  constituted  about  9  percent  of  the  detected 
events  in  that  range. 

E.  REFERENCE  WAVEFORM  FILTER  RESULTS 

Reference  waveform  filters  (RWF)  were  used  on  groups  of 
events  from  three  areas  only:  the  eastern  Kazakh  test  area,  central  Italy 
near  Ancona,  and  a  small  area  midway  between  the  Kamchatka  Peninsula 
and  the  Komandorsky  Islands  .  The  first  two  groups  were  events  which 
originated  from  small  areas  and,  within  each  area,  presumably  had  the 
same  source  mechanism.  The  Kamchatka  events,  except  for  the  reference 
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waveform  event  itself,  were  all  small  and  had  been  undetected  using  both 
bandpass  and  chirp  filters. 


Processing  of  these  events  was  similar  to  the  chirp  filtering 
and  used  the  same  beam  data.  After  applying  the  RWr,  the  output  was  then 
bandpass  filtered.  SNNR  improvements  were  compiled  as  for  the  chirps. 

The  Kamchatka  events  are  listed  in  Table  V-2.  Except  for  the 
reference  event,  none  of  these  events  were  detected  using  the  reference  wave¬ 
form  filter.  The  SNNR  improvement  of  KAM  012*20  when  applied  to  itself  was 
between  4.  5-6.  0  dB  implying  a  magnitude  boost  of  0.  2-0.  3  units.  Allowing  for 
this  gain,  the  probability  of  detection  for  the  m^  =  3.  9-4.  3  range  would  be 
between  40  and  80  percent.*  Included  in  Table  V-2  are  the  results  obtained 
from  ALPA  data  (Strauss,  1973).  ALPA  detected  only  six  out  of  ten 
of  these  events  and  the  surface  wave  magnitudes  of  five  of  these  tend  to  be 
disproportionately  smaller  than  the  given  It  is  likely  that  at  least  a 

few  of  the  events  have  an  erroneously  large  or  are  deep. 

RWF  SNNR  improvements  for  the  Italian  and  Kazakh  events 
are  listed  in  Table  V-3.  The  event  used  as  a  reference  waveform  is  listed 
first.  Kazakh  SNNR  improvements  from  the  RW  F  show  a  distinctly  different 
behavior  from  those  from  the  chirp  filters.  The  vertical  and  radial  Rayleigh 
wave  gains  are  about  2  dB  more  than  ‘he  chirp  gains,  while  Love  wave  improve 
ments  with  the  RWF,  surprisingly,  are  1  dB  less  than  with  the  chirp.  The 
decrease  can  be  attributed  mostly  to  the  event  EKZ*345*04  which  was  actually 
a  double  event  with  origin  times  about  ten  seconds  apart**  and  which  had 
unusually  large  Love  w»vcs.  The  exact  cause  for  negative  improvement  is 
unknown  other  than  that  their  Love  waves  are  obviously  mismatched. 
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TABLE  V-2 


. 


REFERENCE  WAVEFORM  FILTERED  KAMCHATKA  EVENTS 


Event 

m. 

b 

KAM*012*20(R  WF) 

4.8 

KAM* 004  TO 

4.  4 

KOM-009  03 

3.  6 

KAM  '  016  04 

3.8 

KAM  016*11 

3.  9 

KAM  027*20 

3.8 

KAM*032*10 

4.  1 

KAM*  033*  14 

3.  7 

KAM*  05  9*11 

4. 1 

KAM*  05  9*  20 

3.6 

KOM*  0D3*08 

4.1 

TABLE  V- 3 

REFERENCE  WAVEFORM  MATCHED  FILTER  IMPROVEMENT 
IN  SIGNAL+NOISE/ NOISE  RATIO 


PFGIDN  'FK?  '  -  EASTERN  KAZAKH  TF  ST  A  Pf  A 


FVPNT 

r  I  L  TF  R 

RANDPASSFO 

SNNP 

c  NNR 

IMP 

(  DPI 

LOT 

L  P  V 

LPP 

LOT 

LRV 

LRP 

FK  Z / 364/06 

CK73  64/A?04  8 

0.0 

0.0 

0.0 
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0.0 
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FK  7  364 /A  204  8 

- 

6.3 

5.4* 
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The  ten  Italian  events  were  characterized  by  Love  waves  which 
were  virtually  identical,  not  only  in  the  primary  arrival  but  in  the  scattered 
energy  received  up  to  1200  seconds  later  (source-to-array  travel  time  was 
about  450  seconds  and  duration  about  300  seconds'.  This  scattered  energy  had 
decayed  to  background  level  at  1200  seconds  but  was  still  coherent  event-to- 
event.  The  Rayleigh  waves,  however,  were  weak,  confused,  and  relatively 
narrow-band  and  seemed  to  be  of  two  types.  Type  "B"  were  the  waves  from 
IIA  035  09  and  ITA*036*07  while  the  remaining  eight  were  of  type  "A". 

These  two  types  were  almost  identical,  differing  only  by  a  slight  phase  shift  in 
the  waveform  about  midway  through  the  signal.  Both  ’  B’  events  had  exceptionally 
low  (negative)  improvement.  Excluding  the  'B1  events,  RWF  and  chirp  filter 
SNNR  gains  were  about  the  same  for  the  Rayleigh  waves.  RWF  gains  for  the 
Love  wave  were  about  l  dB  larger  than  the  chirp  gains. 

Considering  the  close  similarity  of  the  Love  waves  of  the  ITA 
events,  it  was  puzzling  at  first  as  to  why  the  SNNR  gains  were  not  the  same 
as  for  the  RWF  applied  to  itself.  This  phenomenon  was  also  evidr  it  in  the 
chirp  filtering  results.  It  was  discovered  that  the  method  of  measuring  true 
SNNR,  described  in  part  B  of  this  section,  is  not  always  valid.  The  assumptions 
of  this  method  are  that  the  noise  has  no  signal  energy  present  and  that  the 
bandpass  filter  and  the  RWF  differ  mainly  in  average  filter  amplitude  and 
that  the  difference  can  be  calibrated.  In  many  situations,  these  assumptions 
are  not  true  when  the  noise  spectrum  is  altered  by  the  presence  of  signal 
energy  from  previous  events.  The  non-white  spectrum  of  the  RWF  gives  a 
different  average  output  when  the  spectrum  of  the  input  data  changes  even 
though  the  RMS  level  remains  the  same.  Since  the  filter  spectrum  (in  the 
processing  band)  is  flat  for  the  bandpass  and  chirp  filters  but  not  for  the 
RWF,  the  apparent  gain  of  the  RWF  is  sensitive  to  changes  in  the  noise 
spectrum.  This  sensitivity  results  in  a  RWF  gain  which  is  not  constant  but  is 
dependent  on  the  statistics  of  the  beam  "noise"  at  the  time  of  application.  An 
example  of  this  sensitivity  is  the  behavior  of  the  RWF  on  the  swarm  of  Italian 
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events  beginning  on  4  February  ,c*72.  The  first  event,  used  as  the  RWF,  was 
followed  by  several  large  reported  events  and  probably  also  by  numerous 
unreported  aftershocks.  Table  V-4  gives  the  ratios  of  the  outputs  of  the  RWF 
and  bandpass  filters,  the  chirp  and  bandpass  filters,  aii  the  RWF  and  chirp 
filters  to  both  signals  (Love  wave)  and  noise  for  the  Italian  events. 

The  pertinent  numbers  in  Table  V-4  are  the  means  and  standard 
deviations  (S.  D.  )  of  the  various  ratios.  The  standard  deviation  of  the  ratio  of 
the  RWF  and  bandpass  outputs  on  the  signal  gate  is  small  (equivalent  to  less 
than  0.2  dB).  This  is  true  also  for  the  ratios  of  the  chirp  to  bandpass  and  RWF 
to  chirp.  Conversely,  the  same  ratios  on  the  noise  have  standard  deviations 
three  to  six  times  larger.  Thus  the  variation  in  SNNR  gain  obtained  by  RWF 
(and  chirps  to  a  lesser  extent)  for  essentially  identical  events  is  due  to  subtle 
changes  in  the  noise  spectrum  brought  about  by  remnant  signal  energy,  and 
not  to  any  variation  in  response  to  the  signal  itself. 

The  consequences  of  this  effect  are  that  the  average  gain  of  a 
RWF  should  .lot  be  computed  from  an  earthquake  swarm  and  that  the 
effectiveness  of  a  RWF  may  be  reduced  if  interfering  signal  energy  is  present 
at  the  array  during  the  arrival  of  a  desired  event. 
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TABLE  V-4 

comparison  of  reference  waveform  (ita*035*02),  chirp,  and 

BANDPASS  FILTERS  ON  THE  LOVE  WAVES  OF  THE  ITALIAN* 
EVENTS  COMMENCING  4  FEBRUARY,  1971 


■ 

Rs/  Bs 

Rn/  Bn 

RWF 

Imp 

Cs  /  Bs 

Cn/  Bn 

Chirp 

Imp 

Rs  /Cs 

Rn/Cn 

3  5*02 

1.  904 

1.  274 

1.  495 

1.  172 

0.  944 

l.  241 

1.  625 

1.  350 

35*09 

1.892 

1.  436 

1.  317 

1.  181 

1.  142 

1.034 

1.  602 

1.  257 

35*17 

1.  903 

1.436* 

1.  325 

1.  187 

1. 142* 

1.039 

1.  603 

1.  257* 

3  5*19 

1.  920 

1.  593 

1.  205 

l.  211 

1.  036 

1. 169 

1.  585 

1.  538 

36*05 

1.  909 

- 

1.  213 

1.  075* 

1. 128 

1.  574 

- 

36*07 

1.  932 

1.  573 

1.  228 

1.  187 

1.  075 

1.104 

1.  628 

1.  463 

36*15 

1.  900 

1.  549 

1.  227 

1.193 

0.  951 

1.  254 

1.  593 

1.  629 

3  7*01 

1.  901 

1.  447 

1.  314 

1.  178 

0.  966 

1.220 

1.  614 

1.498 

37*21 

1.  973 

1.  349 

1.  463 

1.  264 

1.  012 

1.  248 

1.  561 

1.  333 

39*12 

1.860 

1.  481 

1.  256 

1.  171 

1.  006 

1. 165 

1.  588 

1.  472 

Mean 

1.  909 

1.  463 

1.  313 

l.  196 

1.  015 

1.179 

1.  597 

1.442 

S.  D. 

0.029 

0.  Ill 

0.  129 

0.  028 

0.  085 

0.  078 

0.  021 

0.  121 

D.  F. 

10 

8 

8 

10 

8 

8 

10 

8 

RMS  noise  level  obtained  from  noise  of  earlier  or  later  event.  This 
data  not  in  statistics 


R  =  RWF  Subscripts: 

C  =  Chirp  Filter  s  Peak  amplitude  in  signal  gate 

B=  Bandpass  Filter  n  RMS  amplitude  in  noise  gate 
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SECTION  VI 

NORSAR  LONG-PERIOD  SUR"  /CE  WAVE  DETECTION  CAPABILITY 

The  detection  history  of  the  383  events  in  the  data  base  was 
used  to  estimate  the  NORSAR  long-period  surface  wave  incremental  detection 
probability  curve.  Smoothed  estimates  of  detection  probabilities  were  ob¬ 
tained  by  a  new  method  (Ringdal.  1974).  This  method  models  detection  as  a 
random  Gaussian  process  and  obtains  maximum-likelihood  estimates  of  the 
process  parameters  based  on  the  experimental  data  at  hand.  For  a  detailed 
description  of  the  theory  and  its  limitations,  and  a  derivation  of  the  likelihood 
functions,  the  reader  is  refered  co  Ringdal  (1974). 

Detections  were  determined  independently  for  each  component 
except  for  the  1971  events  where  detections  were  assessed  only  on  the  vertical 
component.  The  criteria  for  detection  were: 

•  A  peak  in  the  bandpassed  or  matched  filtered  trace  which  was 
3  dB  larger  than  any  other  peak  in  a  15  to  30  minute  time  gate 
centered  at  the  peak  occurence  time.  The  width  of  the  lime 
gate  was  roughly  proportional  to  epicentral  distance  to  account 
for  dispersion. 

•  The  peak  should  occur  within  ±  100  to  +  300  seconds,  depending 
on  distance,  of  the  expected  peak  occurrence  time. 

•  Peaks  should  be  signal-like,  showing  dispersion  appropriate 
for  the  distance  and  region.  For  events  detected  only  by 
matched  filters,  the  matched  filter  output  peak  should  be  pulse¬ 
like  from  the  filter  appropriate  to  the  source  region  and  the 
bandpassed  only  waveform  should  show  some  evidence  of  signal¬ 
like  energy. 
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These  criteria  were  not  absolute  and  for  the  cases  of  threshold 
detection/ non-detection,  judgement  of  the  analyst  played  a  major  part  in  the 
decisions.  In  most  cases,  that  decision  was  conservative  with  the  result  that 
the  false  alarm  rate  was  estimated  to  be  less  than  one  percent  and  that  the  re¬ 
sulting  detection  threshold  estimates  were  also  conservative.  The  detection 
rt  suits  are  presented  in  Figures  V7-1  through  VI-8,  which  include  both  histo¬ 
grams  of  the  event  m^  populations  and  detection  probability  curves.  In  addi¬ 
tion  to  the  observed  detection  percentages  and  the  smooth  detection  probability 
curve,  the  detection  probability  figures  indicate  the  estimated  50  percent  and 
90  percent  detection  m^  levels,  the  90  p  rcent  confidence  limits  of  those  es¬ 
timates,  and  the  standard  deviation  (sigma)  of  the  Gaussian  distribution  fitted 
to  the  observed  data. 

The  detection  behavior  was  observed  also  for  summer  and 
winter  events,  events  in  the  central  Asia  region  and  the  Kurile -Kamchatka  re¬ 
gion,  and  presumed  explosions. 

The  detection  results  by  region  are  as  follows* 

•  The  90  percent  detection  level  for  all  events  combined.  Figure 
VI- 1,  was  at  mb  =  4.  5.  This  is  slightly  better  than  the  esti¬ 
mate  of  m^  *  4.  5  using  just  the  1971  events  last  year. 

•  For  the  Kurile,  Kamchatka,  Okhotsk  region.  Figure  VI-2,  the 
90  percent  threshold  \,as  at  mb  =  4.6.  This  is  slightly  higher 
than  that  obtained  at  ALPA  which  is  45°-55°  closer  than  NORSAR 

•  In  central  Asia,  Figure  VI-3,  NORSAR  obtained  the  90  percent 
level  at  mb  =  4.  5.  These  data  were  must  less  well-behaved  than 
the  Kurile  data  and  thus  the  estimate  is  somewhat  less  reliable. 
In  these  regions,  ALPA  detection  thresholds  are  about  the  same. 

•  The  group  of  presumed  explosions  from  the  eastern  Kazakh  test 
area  was  very  small.  Since  only  two  magnitude  bins  had  more 


VI-2 


DETECTION  PROBABILITY 


MAGNITUDE  C MB  3 


MB  DISTRIBUTION  OF  PROCESSED  EVENTS 
COMBINED  1971  AMD  1972  EVENTS  . 


MAGNITUDE  C MB  3 


MAXIMUM  LIKELIHOOD  DETECTABILITY  CURVE 
COMBINED  1971  AMD  1972  EVENTS 


FIGURE  VI- 1 


NORSAR  LP  SURFACE  WAVE  DETECTION  STATISTICS  FOR  THE 
COMBINED  1971  AND  1972  EVENT  ENSEMBLES 
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FIGURE  VI-2 

NORSAR  LP  SURFACE  WAVE  DETECTION  STATISTICS  FOR  THE 
KURILE,  KAMCHATKA,  AND  OKHOTSK  EVENTS 
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•  FIGURE  VI-3 

NORSAR  LP  SURFACE  WAVE  DETECTION  STATISTICS  FOR  THE 

CENTRAL  ASIAN  EVENTS 
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FIGURE  VI-4 

NORSAR  LP  SURFACE  WAVE  DETECTION  STATISTICS  FOR  THE 
PRESUMED  EXPLOSIONS  FROM  EASTERN  KAZAKH 
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FIGURE  VI-5 

NORSAR  Lr»  SURFACE  WAVE  DETECTION  STATISTICS  FOR  THE 

1972  WINTER  EVENTS 
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FIGURE  VI-6 

NORSAR  LP  SURFACE  WAVE  DETECTION  STATISTICS  FOR  THE 

1972  SUMMER  EVENTS 
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than  one  event,  the  resulting  detection  level  estimates  are 
subject  to  large  error.  With  this  consideration  in  mind,  the 
90  percent  detection  level  from  Figure  VI-4  is  approximately 
at  mfa  *  5.  1  which  is  about  0.4  units  lower  than  for  ALP  A. 

The  1972  events,  after  separation  into  winter  and  summer 
gi cup j,  showed  distinctly  different  detection  behavior  from  season  to  season. 

•  For  the  winter  events,  Figure  VI-5,  the  slope  of  the  incremen¬ 
tal  detection  curve  is  relatively  steep  and  the  90  percent  detec- 

lion  level  occurs  at  m  =4.  5. 

b 

•  The  def  etion  results  of  the  summer  events,  Figure  VI-6,  show 
much  less  consistency  than  the  winter  events.  The  90  percent 

level  is  at  m  =4.2. 
b 

This  seasonal  variation  in  threshold  may  be  significant,  but  is 
partially  caused  by  the  lack  of  adequate  data  for  mfe<  3.  5.  The  slope  of  the 
wintertime  detection  curve  is  rather  steep  indicating  that  only  small  improve¬ 
ments  in  signal-to- noise  ratio  would  be  needed  to  lower  the  90  percent  thres¬ 
hold  bet  that  detectability  of  small  events  drops  very  quickly  with  decreasing 
mb.  Conversely,  the  smalie.  summertime  slope  indicates  that  substantial 
improvements  would  be  reeded  to  lower  the  90  percent  threshold  even  though 
the  probability  of  detection  of  sm,  Uer  events  is  higher. 

The  difference  in  winter  and  summer  detection  probability 
slopes  is  probably  much  less  than  the  indicated  factor  of  two.  The  winter 
statistics  are  relatively  reliable  since  *  ere  are  a  substantial  number  of  events 
at  the  50  percent  level.  However,  for  the  summer  events,  only  a  small  per¬ 
centage  of  the  events  lie  below  the  50  percent  level. 

Matched  filters  are  used  to  increase  the  signal-to-noise  ratio 
of  the  array  output  and  thus  lower  the  detection  threshold.  Figure  VI-7  shows 
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the  population  of  eleven  events  which  were  detected  (on  any  component)  only 

with  the  aid  of  a  chirp  filter.  The  median  magnitude  of  this  group  is  m  =4.0 

b 

which  is  near  the  overall  50  percent  detection  level  of  m.  =  3.  8.  Only  one 

b  1 

event  had  a  magnitude  greater  than  the  overall  90  percent  level  of  m  =  4.  5 

b 

which  implies  that  matched  filters  may  no  substantially  change  the  90  per¬ 
cent  detection  level  but  will  lower  the  50  percent  level.  Eight  of  the  eleven 
were  winter  events  which  would  be  more  positively  affected  by  a  signal-to- 
noise  improvement  by  matched  filters.  Of  the  events  detected  on  only  one 
component,  more  were  detected  using  the  Rayleigh  wave  than  the  Love  wave. 
This  may  be  seen  in  Figure  VI-8.  This  favoring  of  the  Rayleigh  wave  is  due 
in  part  to  the  availability  of  the  radial  component  for  verification  in  a  border¬ 
line  threshold  detection  situation. 
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BEHAVIOR  OF  STANDARD  DISCRIMINANTS 


A.  INTRODUCTION 

The  standard  long-period  surface  wave  discriminants,  M  - 

s 

mb’  AR"mb>  and  AL  -m^,  were  computed  for  each  event  detected.  Measure¬ 
ments  we»e  made  from  the  bandpassed  array  beam  traces. 

The  surface  wave  magnitudes  of  the  Rayleigh  (vertical)  and 
Love  waves  were  computed  using  the  formula: 


or 


where 


M  =  log  (A/T)  +  1.66  log  A  for  A  >  25° 

s 

M  =  log  (A/T)  -f  log  A  +  0.92  for  A  <  25° 

s 

A  =  maximum  peak-to-peak  amplitude  (millicrons)  of 
the  signal 

T  =  period  of  cycle  corresponding  to  A  (seconds) 

A  =  epicentral  distance  from  NORSAR  (degrees) 


The  usual  value  of  1  was  about  22  seconds  with  no  values 
less  than  17  seconds.  There  were,  however,  22  Love  waves  and  10  Rayleigh 
waves  which  had  T  exceeding  25  seconds.  The  majority  of  the  longer  period 
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Love  waves  were  from  Iranian  events  with  periods  of  27-28  seconds.  These 

longer  period  events  were  included  in  the  M  -  m  plots 

s  b 

Nineteen  events  had  either  one  or  both  wave  components 

detected  only  with  a  chirp  filter.  The  M  of  these  events  were  corrected 

s 

using  the  average  chirp  improvement  for  their  area  (Section  V). 

The  Love  wave  magnitudes  were  not  measured  for  the  1971 

events  hence  the  Love  wave  Mg  -  mfe  plots  shown  below  represent  only  1972 
events. 

The  AR  parameter  was  introduced  by  Brune,  Espinosa, 
and  Oliver  (1963)  as  a  measure  of  the  total  Rayleigh  wave  energy  released 
by  a  seismic  event.  The  AL  parameter  is  the  equivalent  measure  of  the 
Love  wave  energy.  The  AR  and  AL  values  were  computed  for  ea  'h  detected 
event  by  summing  the  absolute  values  of  the  beam  outputs  over  a  time  interval 
corresponding  to  a  signal  velocity  range  of  4.  0  to  2.  5  km/sec.  After  scaling 
these  raw  numbers  (Harley,  1971),  the  values  were  normalized  to  a  bodywave 

magnitude  of  5.0  and  a  distance  of  20°.  No  corrections  for  the  ambient  beam 
noise  were  made. 

On  the  following  plots,  earthquakes  are  represented  by 
open  circles  and  presumed  explosions  by  asterisks.  The  values  of  the  various 
parameters,  Mg,  m^,  AR,  and  AL,  are  listtd  in  Table  II-l. 

B-  M  -  m,  DISCRIMINANT 

s  b 

The  Mg-mb  relationships  of  the  Rayleigh  waves  and  Love 
waves  are  shown  in  Figures  VII-1  and  VII-2  respectively.  With  the  exception 
>f  a  few  earthquakes  which  overlap  into  the  presumed  explosion  group  there 
is  good  separation  between  earthquakes  and  presumed  explosions  for  both  Love 
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FIGURE  VII- 1 

M  VS  m  DISCRIMINANT 
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COMBINED  1971  AND  1972  EVENTS 


VII-3 


1 


3  4  5  6  7 


MB 

FIGURE  VU-2 

M  VS.  m.  DISCFIMINANT  1972  EVENTS 

S  D 


VII -4 


and  Rayleigh  waves.  The  overlapping  events,  indicated  by  tags,  will  be  dis¬ 
cussed  below.  Minimum  separation  between  classes  is  about  0.4  M  units 
for  Rayleigh  waves  and  may  be  slightly  more  for  Love  waves.  There  is  a 
distinct  increase  in  the  slopes  of  the  Mg  trend  lines  of  the  earthquake  points 
at  mb>  5.  Similar  curvature  cannot  be  seen  for  the  presumed  explosions 
because  of  insufficient  data.  A  lack  of  both  large  magnitude  earthquakes  and 
low  magnitude  presumed  explosions  prevents  conclusions  as  to  class  separa¬ 
tion  at  the  low  and  high  magnitude  ranges.  In  general,  separation  is  not  as 
good  as  reported  last  year. 

The  events  tagged  in  Figures  VII- 1  and  VII-2  are  listed  below. 

TABLE  VII- 1 

EVENTS  HAVING  POOR  SEPARATION  WITH  M  -  m  DISCRIMINANT 

s  b 


Tag 

Name 

m, 

b 

NORSAR 

ALPA 

M 

s 

M 

s 

M  * 
s 

M  * 
s 

Rayleigh 

Love 

Rayleigh 

Lo  ve 

i 

IWP*028*10 

3  9 

3.  9 

—  i 

2 

KAM*1 99*08 

B 

3.  2 

4.  2 

3.9 

4.  3 

3 

GRE*190*05 

B 

2.  7 

2.  7 

ND 

ND 

4 

YUG*180*01 

4.  6 

2.  6 

3.  1 

ND 

ND 

5 

BSA/210/19 

4.  5 

2.  6 

— 

— 

6 

KUR*077*07 

5.  2 

3.  5 

3.  1 

3.  8 

3.9 

7 

IRA  *18  5*0  2 

5.0 

4.  0 

2.  8 

— 

— 

8 

I TA  *0  35*04 

4.  8 

ND 

2.  1 

ND 

ND 

*  Magnitudes  obtained  at  ALPA.  Special  Report  No.  8 
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These  events  overlap  into  the  population  of  presumed  explosions 
for  the  following  reasons: 


•  Events  1  and  8  had  values  assigned  by  only  one  station 
each  of  which  were  at  non-teleseismic  distances  (A<  15°). 

•  Events  2  and  7  each  had  only  one  component  with  a  low  M 

s 

probably  due  to  radiation  pattern  effects. 

•  Events  3  and  4  were  from  the  same  area  and  had  many 

stations  reporting  These  may  be  deep. 

•  Event  5  was  reported  only  by  LASA  and  was  near  LASA's 
range  limit.  The  m^  is  probably  erroneous. 

•  Event  6  had  values  measured  at  ALPA  which  were  not 

unusually  low.  We  have  no  explanation  for  the  low  M  's 

s 

measured  at  NORSAR. 

Ms  *  mb  was  plotted  on  a  regional  basis  for  the  Kurile- 
Kamchatka  area.  (Figures  VII-3  and  VII-4)  and  the  Central  Asian  region 
(Figures  VII-5  anu  VII-6).  The  presumed  explosion  ensemble  is  included 
for  comparison. 

The  Kurile-Kamchatka  -  m^  data  appears  to  show  good 

separation,  except  for  the  tagged  events,  for  both  the  Love  and  Rayleigh  waves. 

The  range  of  mb's  of  the  earthquakes  and  presumed  explosions  barely  overlap 

each  other  thus  true  separation  for  m,  >  5.  2  is  not  known. 

b 

The  Mg  -  mfe  data  for  Central  Asia  show  distinct  differences 
between  Love  and  Rayleigh  modes.  Separation  is  good  in  both  cases  and  a- 
mounts  to  about  one  magnitude  unit  for  the  Rayleigh  data.  The  somewhat  greater 
scatter  of  the  Love  M^'s  tends  to  reduce  the  distinctness  of  separation  although 
separation  is  still  good. 
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FIGURE  VII-6 

Mg  VS.  mb  DISCRIMINANT  REGIONS 
PAK,  TKZ,  TIB,  CHI,  SIN,  EKZ 
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(  .  AR-mb<  AL-mb  DISCRIMINANTS 

The  plots  of  AR  and  AL  are  shown  in  Figures  VII-7  and  VII-8. 
With  the  exception  of  one  event  with  a  poor  ny  IWP*028*10,  there  is  good 
separation  between  earthquakes  and  presumed  explosions  with  both  AR  and 
AL.  Minimum  separation  is  by  a  factor  of  four.  The  largest  AR  from  a 
presumed  explosion  was  produced  by  WKZ/356/06  with  that  for  EKZ*345*04 
almost  as  big.  The  latter  was  actually  a  double  event  with  origin  times  ap¬ 
proximately  eight  seconds  apart  and  produced  the  largest  AL  value  of  its 
group.  This  AL  was  three  times  larger  than  the  next  biggest  AL  in  its  group 
and  cannot  be  explained  simply  on  the  basis  of  constructive  interference  of 
two  otherwise  normally. sized  signals. 

AR  and  AL  show  no  significant  advantages  of  one  over  the  other 
even  though  AL  might  be  expected  to  show  slightly  better  separation  because 
of  reduced  excitation  of  Love  waves  by  presumed  explosions.  Rayleigh  wave 
energy  from  presumed  explosions  apparently  is  also  relatively  small. 

AR  and  AL  are  closely  related  to  Mg  because  each  is  a  mea¬ 
sure  of  the  seismic  energy  released  by  an  event.  Thus  AR  and  AL  measure 
area  rather  than  a  peak  amplitude  and  period,  they  would  tend  to  be  less  in¬ 
fluenced  by  small  quirks  in  the  waveform  dispersion.  They  also  appear  to  be 

less  sensitive  then  M  to  errors  in  m,  values. 

5  b 
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SECTION  VIII 


CONCLUSIONS 


A.  MAJOR  RESULTS 

The  major  conclusions  from  the  five  studies  of  the  NORSAR 
evaluation  are  summarized  below.  These  conclusions  are  based  on  the  overall 
results  obtained  (luring  the  past  two  years. 

1 .  Noise  Anal ysis 

•  The  ambient  noise  at  NORSAR  is  strongly  seasonal  with  summer¬ 
time  levels  of  5-10  m/x  RMS  (20-40  seconds)  and  wintertime 
levels  of  8-15  m^x  .  The  noise  appears  to  be  almost  entirely 
made  up  of  propagating  surface  wave  energy. 

•  During  the  winter,  noise  levels  often  rise  to  levels  greater  than 
25  m fi  RMS.  These  increased  levels  are  definitely  caused  by 

st  '’■ms  in  the  North  Atlantic  Ocean,  however,  the  noise  genera¬ 
tion  mechanism  is  not  understood.  The  duration  of  these  high 
levels  may  he  between  1 2  to  36  hours. 

•  The  source  directions  of  the  ambient  noise  are  also  seasonally 
dependent  and  seem  to  have  a  few  preferred  azimuths.  In  the 
winter,  noise  is  primarily  from  the  north  and  the  west.  In  the 
summer,  the  primary  noise  azimuth  is  east.  "Summer"  noise 
can  V  differentiati  d  from  "winter"  noise  by  the  change  in  noise 
azimuth  from  west  to  east  and  vice  versa.  The  direction  changes 
are  not  strongly  correlated  with  noise  levels. 
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2. 


Array  Processing  Performance 


These  conclusions  are  based  on  results  from  eleven  winter 

noise  samples. 

•  Within  the  MCF  design  gate,  the  MCF  processor  achieved  2  to 

8  dB  more  reduction  in  noise  than  the  BS  processor  in  the  signal 
processing  band  of  0.025-0.  059  Hz.  Outside  of  the  design  gate, 
this  advantage  dropped  slightly  to  0  to  6  dB.  Best  MCF  gains 
were  not  obtained  when  noise  levels  were  highest. 

•  Array  gains  of  both  the  MCF  and  BS  processors  were  seriously 
degraded  when  array  beam  look  directions  were  close  to  the 
predominant  noise  azimuths.  For  look  directions  significantly 
different  from  the  noise  directions,  the  array  gains  improved 
and  the  MCF  obtained  2-5  dB  additional  gain  over  the  BS  for 
on-design  noise  and  0-4  dB  for  off-design  noise. 

•  The  MCF  typically  achieved  3  dB  more  signal -to-noise  ratio 

improvement  than  the  BS  processor  in  actual  signal  extraction. 

This  suggests  than  an  MCF  processor  could  lower  the  NORSAR 

detection  threshold  in  winter  by  0.  15  M  units 

s 

•  In  a  comparison  of  MCF  and  BS  processors  using  the  full  array 
and  a  smaller  subarray  configuration  ,  the  MCF  over  BS  array 
gain  improvements  were  slightly  larger  for  the  subarray. 

3.  Matched  Filtering  Performance 

•  A  matched  filter  parameter  was  highly  successful  as  the  basis 
for  defining  seismic  regions  in  Eurasia.  The  chirp  filter  length 
giving  optimum  signal -to -noi se  ratio  improvement  was  found  to 
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4. 


show  regular  and  consistent  variations  for  various  regions. 
Maps  were  produced  showing  contours  of  equal  chirp  length. 

Average  chirp  filter  SNNR  gains  for  different  regions  ranged 
from  0.  7  to  4.  8  dB  with  typical  gains  of  2  dB.  The  variation  in 
gain  within  a  region  was  rather  large.  The  standard  deviation 
of  gain  for  the  regions  ranged  from  0.  5  to  2  times  the  average 
gain  with  usual  values  of  0.  6  of  the  gain. 


Chirp  filters  produced  a  small  increase  in  the  number  of 
detected  events.  The  m^'s  of  these  events  were  in  the  neigh¬ 
borhood  of  the  50%  detection  level.  These  additional  detections 
lid  not  significantly  affect  the  90%  detection  threshold. 


The  SNNR  gain  of  refeience  waveform  filters  was  shown  to  be 
degraded  in  the  presence  of  signal-like  energy  in  the  background 
noise.  This  implies  that  a  RWF  may  not  be  effective  in  separating 

a  desired  event  from  an  interfering  event.  Chirp  filters  are  less 
adversely  affected. 

NORSAR  Surface  Wave  Detection  Capability 

All  detections  were  measured  with  a  false  alarm  level  of  less 

than  one  percent.  The  90  percent  detection  thresholds,  by  region, 
are  the  following: 

All  events  combined:  m  =  4.  5 


Kurile,  Kamchatka,  Okhotsk:  m  =  4  6 

b 

Central  Asia:  m,  =  4.  5 
b 

Eastern  Kazakh  (presumed  explosions  only):  m  =5.1 

b 

The  90  percent  detection  threshold  at  NORSAR  shows  a  strong 
seasonal  dependence  .  This  may  be  attributed  to  the  increased 
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ambient  noise  level  during  the  winter.  The  winter  threshold  is 

at  m  -  4.  5.  The  summer  threshold  is  at  m  =  4.  2. 
b  b 

•  Chirp  matched  filters  are  probably  not  very  effective  in  lovering 
the  90  percent  detection  threshold. 

5.  Behavior  of  Standard  Surface  Wave  Discriminants 

•  Reasonably  good  separation  was  obtained  by  the  M  -  m 

s  b 

discriminant.  Overlap  of  event  classes  was  found  for  eight 

events,  however  none  were  overlapped  on  both  Rayleigh  and 

Love  wave  M  .  There  was  essentially  no  difference  in  sepa- 
s 

ration  between  Rayleigh  M  -  m,  and  Love  M  -  m  . 

s  b  s  b 

•  AR  and  AL  appear  to  be  better  discriminants  than  M  at 

s 

NORSAR.  Good  separation  was  obtained  for  both  AR  and  AL. 
With  the  exceptioi  of  one  presumed  explosion  (actually  two 
events  with  origin  times  eight  seconds  apart),  minimum  sepa¬ 
ration  by  AL  was  by  a  factor  of  four. 

B.  PLANS  FOR  FUTURE  ANALYSIS 

The  results  from  the  NORSAR  evaluation  to  date  have  suggested 
studies  in  the  following  areas: 

•  Continue  to  monitor  the  ambient  noise  field  at  NORSAR  with 
particular  emphasis  on  seasonal  changes. 

•  Investigate  further  the  array  processing  gains  of  MCF  and 
beamsteer  processors  using  both  the  full  array  and  a  subarray 
for  both  signals  and  noise. 
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Refine  the  boundaries  of  seismic  regions  defined  by  matched 
filter  behavior  by  expanding  the  event  data  base. 

Use  the  expanded  data  base  to  improve  estimates  of  the 
NORSAR  detection  threshold  and  discrimination  capability. 
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